SOLUTION OF LAMINAR COMPRESSIBLE BOUNDARY 
LAYER EQUATION BY PARAMETRIC DIFFERENTIATOR 

METHOD 


BY 


AJOY KUMAR MITRA 


He 

H 


DEPARTMENT OF MECHANICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

FEBRUARY, 1970 



CENTRAL LIBRARY 
Indian Institute of Technology 
KANPUR 
yWesyCi 

Class No. . 6 . 2 . 3 .-. 1 . 331 ... 
ri£ 97-4 



SOLUTION OF LAMINAR COMPRESSIBLE BOUNDARY 
LAYER EQUATION BY PARAMETRIC DIFFERENTIATION 

METHOD 


A Thesis Submitted 
In Fulfilment of the Requirement 

for the Degree of 

MASTER OF TECHNOLOGY 



BY 

A JOY KOMAR MITRA 



SOL. I’- 

f-i an 

DEPARTMENT OF MECHANICAL ENGINEERING 



The author recognises the invaluable suggestions 
and the precious help of Dr. Oberai in selection 
and completion of the problem. 




CERTIFICATE 


This is to certify that the present work has been carried 
out under my supervision and has not been submitted elsewhere 
for a degree* 



/r ). 

/ 

( Madan Mohan Oberai) 

Assistant Professor 
Department of Mechanical Engineering 
Indian Institute of Technology , 
Kanpur, IHDIA. 


\ 

1 


This thesis has been approve 

Master ui »■ ‘‘‘W 1 

i„ accoteU.ee with the 

regulations of the Indian 
Institute of Technology Kanpur 
TVited. C.^' ?• 



ABSTRACT 


In the present work, two -dimensional laminar com- 
pressible boundary layer problem under the assumption of 
either Prandtl number equal to unity or that of low mach 
number, has been solved numerically using Parametric Diffe- 
rentiation Method. At the end, a general program based on 
Fortran IV has been developed by which through the use of 
the initial value obtained by the Parametric Differentiation 
Method as starting guess, the original non-linear equation 
can be solved to any degree of accuracy. 
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CHAPTER 1 


INTRODUCTION 

The equations governing the two dimensional compres- 
sible laminar boundary layer are coupled in addition to being 
non-linear. Except for the special case of a flow with Prandtl 
number unity over an insulated surface, these equations are 
difficult to solve. 

With no standard method can the solutions of these 
two point boundary value problems be expressed in a closed 
form. The following two methods are generally resorted to 
for solving these equations. 

1. Forward integration 

2, Integration by successive approximation 

Both of these methods are highly iterative and a good 
amount of guess work has to- be done for rapid convergence to 
the solution. The aim of the present work is to dispense 
with the guess work and to solve the equations by parametric 
differentiation (Ref. 5). As we shall see, this method reduces 
the equation to b,§, solved to a linear system which, though of 
higher order, can be solved by a predetermined number of 
iterations. 

Satisfactory results obtained in the case of the com- 
pressible boundary layer encouraged us to develop^ a general 
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program which can handle all those ordinary differential 
equations which can he solved hy parametric differentiation 
method. ' In the last chapter we have presented this general 
program. Due to the flexibilities incorporated in the pro- 
gram, it can handle, in addition to the reduced linear sys- 
tem of equation, the parent non-linear system of equations. 
This flexibility is used when a high' accuracy is needed. Here 
the solution is obtained in two steps in the first step, the 
usual solution is obtained using parametric differentiation 
method. In the second step, the nonlinear equations are 
directly solved. Here the missing initial boundary conditi- 
ons are taken from the solutions obtained in the first step. 
As initial boundary condition taken are very close to the two 
one. This method converges very rapidly in a few iterations 
to the ti.ue solution . 



CHAPTER 2 




PORMULATION OP THE PROBLEM (Ref. 4) 


2.1 Stewart son’ s Equation: 

i 

Steady two-dimensional compressible laminar boundary 
layer for .perfect gr - "• are given by 


~5F ( pu ) + Ty 


pu li + ^ n 


0u 

y 


H + ^ If 


( Pv) = 0 (Cont. Eqn. ) 

: ~ + Jy 0 U §^) (Mom. Eqn. ) 

9p 


(D 

( 2 ) 


u ^ + JL ( A 

9x 3y ^ Pr 9y 


) + 


a r 

(Energy Eqn. ) (3) 


where , 

p = pressure 
P = density 

u = longitudinal velocity 
v = transverse velocity 
x = longitudinal coordinate 
y = transverse coordinate 
h . = enthalpy 

h ~ coefficient of viscosity 

Pr = Prandtl Humber 

The viscosity law to be assumed is 



( 4 ) 
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The constant /\ is used to match the viscosity with the 
Sutherland value (Ref. \3) at a desired station. If this station 
is taken to be the surface assumed to be at constant temperature 
the result is 

a = rrrTTT ( ) (5) 

v ' w o' v t + K ' 
w su . 

where 

t = wall temperature 

t Q = free strain stagnation temperature 
K = Sutherland's constant 

SIX 

2.2 Stewart son' s Transformation: 

Velocities in the equation of motion (1) to (3) can be 
expressed in terms of derivatives of stream function as 


f y = ( PV/ f 0 ) 

Vi = ( fo 5 


Now introducing transformation 

a 


X 


X P Q 
r _e 

j p 

o *0 


dx 


• 3s 


fo 


dy 



( 7 ) ' • 


(Subscript e refers to local condition at the outer edge of 
the boundary layer (external). The subscript o refers to the 
free-stream stagnation values.) Where a is sonic velocity. 



J 


We get the following equations: 


Uy + Vy = 0 


UU 


I 


VU V = u U Y S + ^ tt 


e ex 
S 


o YY 


US X+ VS Y = (-g 


1-Pr 
Pr 1 


((y-i)/2)a 


f 2 ■, 

e / 


1 + ((y -1 )/2 )M 


U e 


Y Y 


( 8 ) 


(9) 


where 

Y = ratio of specific heats 

-h = * 

a = h s /h Q (10) 

a is local sonic velocity. 

s is enthalpy function 

and h is local stagnation enthalpy, 
s 

Transformed velocities U, V are related to stream func 
tion hy the following equations: 

U = -Y y , V = if. x 

The transformed longitudinal velocity U is related to the 
longitudinal velocity in physical plane by 
H = (a Q /a e ) u 

The boundary conditions applicable to eqn. (7) to (9) are 
TJ(X, 0) = 0 

V(X, 0) = o 

s(X ? 0) = s w or ^ (X, 0) = (||)^ (11) 

lim s = 0 

Y-oo 
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Now under the transformation 


f = f(T) ) ri'2 4 0 U e X)/(m+1) 
r ! = Y fo: 5 (m+17'U e 7ri^xT 


( 12 ) 


The above system of partial equations is reduced to following 
ordinary diffrential equation 


f +ff = p(f - S) 


S + Pr f S = ( 1-Pr) 


(V- 1)M? 


| 1 + ( (4-1 )/2)Mt 


, 1 tft HO 

(f f + f • ) (1 


The pressure gradient parameter |3 is defined as 


P = 

where m is given by 


2m 

m+1 


u e = Cx“ 

where C is- constant. 

The velocity ratio U/U = u/u = f* 

In above equation prime denotes differentiation with 
respect to . 

The boundary conditions are 

f(0) = f'(0) = 0 

S(0) = s 

v 1 w 

t 

Lim f = 1 


( 14 ) 
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Now in the energy equation, the right handside is not 
functionally consistent with arbitrary Mach no. T r e which is 
a. function of X. 'Hence to be consistent with left handside 
right handside - of energy equation must be either constant 
or function of¥| . Now this is true for following cases 

1 . ’Then external Mach number is constant as is in the 
case of flate plate. 

2. External Mach number is negligibly small. 

3. Prandtl number is equal to unity 

4. Specific heat ratio is equal to unity. In reality !■,. 
cannot be equal to uhity. 

5. Mach number is -very high (Me - °°) as in this case 
right hand side will be equal to 2. 

Here problem will be treated assuming either Prandtl 
number is equal to 1 or Mach number is very small. 

Hence system of equations to be solved are 

Hi " M . » O v , . ^ x 

f+ff=p(f-s) (15) 

s + fs Pr = 0 • ( 16 ) 

with boundary conditions (14) 



CHAPTER 3 


METHOD OE SOLUTION 

3. 1 Reduction of Nonlinear Equations to Linear Form: 

Essence of method "based on parametric differentiation 
'is to differentiate the original equation with respect to a 
parameter, and to solve the resulting linear equation. 

Now to apply this method we have to know the solution 
of the equation for one value of parameter and then we inarch 
forward for other values of the parameter. 

Now differentiating equations (15) and ( 1 6 ) with res- 
pect to (3 we get, 

G + Gf + fG = (f ^-S) + p(2f G -T) (17) 

t" + Pr ^(fT 1 + GS') = 0 (18) 

with boundary conditions 

G*(0) = G(0) = 0, T(0) = 0 , T (°°) = 0 

where 

df/dB = G 

(19) 

dS/dp = T 

where prime denotes differentiation with respect to 

Equations (17) and (18) are linear in G and T. 

Now after solving for f and S for initial value of 
p, we substitute, for f and S in linear equation (17) and (18) 
and evaluate G and T. From these values of G and T, we 



9 


solve for f and s from equations (19). This cycle is repeated 
till (3 reaches the desired value. Both the systems of equa- 
tions (17, 18) and (19) are solved using Runge-Kutt a method 
with proper increment in Yj and (3 respectively. 

Initially such a value of parameter is chosen that at 
this value either the solution is known or it can he obtained 
very easily. 

Here the initial value of parameter p is chosen as zero. 
For this value of p, the equations (l 5, 16) are reduced to 

f + ff = 0 (20) 

s” + Infs' =0 (21) 

with boundary conditions (14). 

Equation (20) is Blasius equation whose solution is 
known .or can he obtained very easily. 

Solution for other values of (3 from that at the initial 
value of |3 is obtained following the method outlined above. 

In actual numerical solution, increment in 1] has been 
taken equal to .05 and that in p. equal to 0.1 for favourable 
pressure gradient and -0.05 for adverse pressure gradient. 
Solutions were obtained for various Prandtl numbers and wall tem- 
peratures at different values of p. 

s We have seen equations to be solved for p equal to zero 
are uncoupled and simple but still it is nonlinear. We can 

"Pn r*"h"h t* Qinrml i -Ptr irm'hl om "htr ■fnY'+ViP ’n ncjo n*F ‘nflrflmp'h'ri (7 
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differentiation method. In Blasius equation there is no 
parameter hut we can artificially introduce a parameter as 
shown helow. 

Let us consider the equation, 

f"’ + (l-P)f' + P ff" = 0 (22) 

with boundary condition 

f(0) = f' (0) = 0 

f (°°) = 1 

Above equation for P equal to one is reduced to Blasius 
equation. For P equal to zero equation is reduced to, 

H ! ?! 

f + f = 0 (23) 

with boundary condition 

f(0) = f' (0) = 0 f* (oo) = i. 

Solution of equation (23) is very simple and is given by 
f = ^ + e -V "l - 1 

Now differentiating equation (22) with respect to P, 
we get, 

w"'+ ( 1-P) w" + f" + ff” + PWf” + Pfw" = 0 .(24) 

with boundary conditions 

W(0) = w'(o) = 0 

W(°°) = ■ 1 


( 25 ) 
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Vlhe re , 

df/dP = W (26) 

Equations (22, 24 and 26) are solved in the same way as 
equations (17, 18 and 19). 

Equation (23) has been solved with increment in 
equal to 0,05 and that in P equal to 0.05. Results obtained 
tally fairly with the standard result. Results are shown in 
table (1) oar the next' page. 

Similarly initial solution for both the equations 
(15, 16 ) can be obtained by assuming suitable equations. 

This has been illustrated in ^Chapter 5. 
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Table 1 

Solution of Equation (22) by Parametric Differentiation 
for Pr = 1 , (Numerical solution for Blasius Equation) 


ETA (V]) 

P 

ETA (f |) 

F 

0 

0 

3.8 

2.583 

0-2 

.0094 

4.0 

2.782 

0.4 

.0375 

4.2 

2.982 

0.6 

.0843 

4.4 

3.181 

00 

• 

O 

.1495 

4.6 

3.381 

1.0 

.2328 

4.8 

3.581 

1.2 

.3333 

5.0 

3.781 

1.4 

.4503 

5.2 

3.981 

1.6 

.5824 

5.4 

4.181 

1.8 

. 7282 

5.6 

4.380 

2.0 

.8860 

5.8 

4.580 

2*2 

1.054 

6 . 0 

4.780 

2.4 

1.231 

6.2 

4.980 

2.6 

1.414 

6.4 

5. 180 

2.8 

1.602 

6 . 6 

5.380 

3.0 

1.794 

6 . 8 

5.580 

3.2 

1.989 

7.0 

5.780 

3.4 

2. 186 

7.2 

5.980 

3.6 

2.384 

7.4 

6.178 



CHAPTER 4 


RESULTS OBTAINED 

4.1 Velocity and 'Enthalpy Function: 

The velocity and enthalpy functions are presented in 
tabular form. 


The distance y normal to the surface in the physical 

plane is related to the similarity variable Tj through 

equation (6) and (12) as (Ref. 4). 

p • a 
1 o _c 

' a 


y 


, I 
Vo 


e -4 m+ 1 


Ti-t 

TT f t " d,! - 

e o o 


where 


t/t = (1 + 


‘2-1 


M 2 ) S - M 2 f' 2 
e J 2 e 


(27) 


2 e * 2 

4 . 2 lntegral Thicknesses: 

The boundary layer integral thicknesses in the 
transformed plane are defined by the following relations (Ref. 4 

Displacement thickness: 


tr 


X 


0 


tr 

x 


Thermal thickness : 


\ 133+ 1 

U X 
e 


s| 

O 1 

ckness: 


I 

1 m+1 

V 

4 2 

^0 


/ (S - f ) <*T) 
0 


(28) 


/ f (1 - f ) drj 
o 


(29) 


i m V 
2 4 W 


/ s dT| 


(30) 
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Convection thickness: 



m +1 

2 ” 


V 

N 0 


I (s-i ) f’d = -s' 

o 


(31) 


Numerical values of these thicknesses a-re evaluated for 

different values of 8 , Pr and S . 

w 

4.3 Shear and skin friction: 


The quantity that is of primary interest in boundary 
layer calculation is the shear stress at wall '~C W which in 
nondimensional form can be given as (Ref. 4) 



1 

2 


'Vw 

fw u e 



4' (2 Mi+s w )) 


\ mjj j) o 
1 



which can be rewritten as 


C iTRe 


w _ 


w 


m+i d I n X 
2 "d In" x 


(33) 


where is called local skin friction coefficient. 


Re 


V 


w 




w 


w subscript denotes properties at wall. 


4.4 Heat Transfer: 

S* = ds/dv' at wall corresponds to heat, transfer 
across boundary layer. This is related to stagnation enthalpy 
derivative in the physical plane by 

\ ro +1 


h. 


(r-s) 


( 


f a 


e 

Tv 


) s 


(34) 


15 


A non-dimensional quantity in connection with heat 
transfer can he introduced (Ref. 4) as 


Ntt = 


x( 6 t/ 6 y) 


w 


( ~w N irs - ” \ m +1 d In X 
g ^-_ ^ v Ke w 4 I ” d In x 


t - t 
o w 


(35 


Reynolds analogy: 


A simple modified Reynold analogy can he defined as 


^f Re w 

“ThT"' 


(-s. 


2 fw 

7 & 


(36 


w 


• 1 ) 


This quantity is the reciprocal of usual Reynold analogy 
quantity. This factor is tabulated for different values of 
p at different values of Pr and Sw . 

Numerical results are given in Appendix 1 . Computer 
program for this is given in Appendix 2. 


4.4 Some Comments on the Result; 

In Pig. (1, 2, 3)velocitv profiles have "been shown. We 
observe that 

i) for a given wall temperature and Prandtl number «, the 
initial slope of velocity increases with the pressure gra- 
dient parameter (Pig. 7 also). 

ii) for adverse pressure gradient, an inflexion occurs 
within the boundary layer that moves outwards with increase 
in magnitude of negative pressure gradient parameter. 
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iii) for Sw > 1 , (or for heated surface) under favourable 
pressure gradient a velocity overshoot (Fig. 2) occurs within 
the boundary layer. This overshoot increases as the pressure 
gradient parameter becames more favourable. The over shoot 
takes place due to the fact that when wall is heated in a 
favourable pressure gradient, the density in certain region 
of boundary layer is lowered so that in spite of viscous drag, 
the flow there is accelerated more than the external flow by 
the pressure force. 

iv) velocity boundary layer increases with increase in 
Prandtl number. 

v) In Fig. (3 , 4 , 5) stagnation temperature profiles have 
been drawn. We see ' for a given Prandtl no. and specified 
wall temperature profile vary very slowly for different pres- 
sure gradient parameters This variation becomes prominent 
with increase in the wall temperature. Like initial velocity 
gradient , slope of initial temperature gradient at wall in- 
creases as the pressure gradient becomes more favourable 
(Fig. 8 also). 

vi) Thermal boundary layer thickness decreases with in- 
crease in Prandtl number. 

By inspecting tabulated results, we observe that 
for low wall temperature in favourable pressure gradient 
displacement thickness is negative. This is due to the fact 
that the fluid in contact with the cold wall has higher den- 
sity than that in the external flow. Hence more fluid /area. 



CHAPTER 5 


SOME RETAILS OP COMPUTER PROGRAMMING 

5-1 Iteration Scheme Used for Satisfying Boundary Conditions 

Before explaining about the iteration scheme, we will 
mention one of the important properties of the ordinary linear 
differential equation. 

Let the system of linear differential equation be given 
by 

1(D) = 0 

How this can be reduced to a set of linear differential equa- 
tion of order one, number of equations of order one being equal 
to the suth of the orders of the original system-- e-f equations. 
Any derivative of any dependent variable in the original equa- 
tions will correspond to one dependent variable in this new 
set cf equations of order one. 

Let x 2 , Xj, ... be dependent variables in the new 
set of equations. Let x^(0) , x^(0),... x (0) are missing ini- 

tial conditions whose value we have to know for solving the 
problem as initial value problem. , 

Let the conditions to be satisfied at the end point be 

x a (R) = M(1) 

x c (R) = M(2) 
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where FN(l), M(2), ... EN(n) are values of variables x , x etc, 

<3, C 

to be satisfied at the end point R. Let us relabel the missing 
initial conditions as G-( 1 ) , G(2), . .. .G(n). 

Thus 

G(1) = x.(0) 

G(2) = x^(O) 


G(n) = x p (0) 

Similarly, x , x , ... x are : yelebclled as 

d O , j 

Y(1) = x a (R) - M(1) 

Y(2) = x 0 (R) = m(2) 

Y(n) = x q (R) = M(n) 

Let G(i) = C(i), for i = 1, n 

where C(i) is a value assumed for G(i). 

Now for linear equation it is well known that if all the 
initial conditions are kept fixed, and an increment is given to 
the initial value of a variable then 

AT.(j) = 0. ± &G(i) (57) 

where AT^(j) is increment in Y(j) due to the increment G(i) 
in G(i), all other initial values of variables being kept constant 
0.. is a constant. 

Graphically it is shown in the adjacent figure 
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Now if there are increments x^ , x^ in G(l), G(2) G(n) 

then total increment in T(j) will he given by 

AY(d) = S AT, (3) = S c x. (38) 

i=1 1 i=1 31 1 

Now we can treat the problem inversely in the following manner. 
Giving different values to initial conditions, we can calculate 
coefficients like C,, from (37), 

Now if A7(o) is departure of Y(j) from the prescribed IN(3), 
when initial values are assumed as 

G(i) = C( i) , i = 1 , n 

Now for a given AY( t l), we can solve for x^ from the 
simultaneous algebraic equations (38). After obtaining x^, correct- 
initial conditions which will satisfy the given boundary condition 
will be given by 

G(i) = C(i) +■ x(i)y i = 1, n (39) 

Now above is strictly true for linear equations but 
approximately true’ for non-linear equations where guessed initial 
values are near the true values. 

In above method of interpolation for obtaining correct 
initial values, effect of modifications of all initial values 
on any particular final value has been taken into consideration. 
Thus in non-linear case it can be hoped that convergence to the 
correct' initial value will be obtained more rapidly in the scheme 
of , interpolation than that can be obtained in a scheme where .into 
polation is done only between one initial value and corresponding 
■final val up. 
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5 . 2 Example 

General program can best be understood with reference to 
some equations to be solved. Let us take the equations (15) 
and (16) as an example. 


( 40 ) 


Consider the following equation 
f'" + B (ff" + (s - f' 2 ) p) + f” ( 1-B) = 0 
s M + B fs' Pr + s' - Br* = 0 
with boundary conditions f ' (0) = f(0) = 0, f * (°°)= 1 

s (0) = Sw, s(°°) = 1 

Above set of equations is reducible to equations (15? 16) 

for B ^equal to unity. Now for B = 0, equation (40) is reduced 

to a very simple form (41) 

ftn + f M ( 1 -B ) = 0 

s " + s ' =0 

whose solution is given by 

f = v* + e~^ -1 1 

V. 

s = 1 + (Sw - 1 ) e~ ^ 

Erom the solution of equations (40) for B equal to ^zero? we 


(41) 


(41a) 


can march forward by parametric differentiation method to get 
the solution at B equal to unity. 

Though equation (40) looks more complicated than the 
parent equation nonetheless it is chosen as it gives the start in 
solution readily. 

: Now differentiating (40) with respect to B we get 
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G ,M +(ff , +(s-f' 2 ) +B(Gf”+fG"+(T-2f' )p)-f"_G"( 1-B) =0 


BT'-s' = 0 


(42) 


(43) 


1"+ fs' Pr + B Gs' Pr + BfT'Pr + T' - 
with boundary conditions 

G(0) = G ' (0) = G'(~) = 0 
and T(0) = T(°°) = 0 

where 

r _ df m _ ds 

dB ’ 1 ~ dB 

Parent equations are 

f m + ff" + (s-f ,2 )p = 0 

s " + fs ' Pr = 0 

with boundary conditions 

f(0) = f ' (0) = 0 
f'(°°)= 1 

and s(0) = Sw, s(°°) = 1 

Now we have 3 sets of equations. 1st set is equation (41), 

second set is (42) and third one is eqn. (43). In PDM (Parametric 

Differentiation Method), we first solve 1st set, then second set. 
If higher accuracy is required we solve the third set of equatiom 
directly. In solving the third set we obtained , starting values 
by solving the second set of equations by PDM. 

In this program any dependent variable is represented as 

• f 

double subscripted variable Q(n, m) where n denotes value of 
the independent variable at which dependent variable Q is beihg 
considered, and m denotes different dependent variables. Corresnon- 
dence between m and different dependent variable is done in a 
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systematic way. With the order in which equations are written 

in & 2) program assigns 

G = Q(n, 1) 

G' = Q(n, 2) 

G" = Q(n, 3) 

G’" = Q(n, 4-) 
and I = Q(n, 5) 

T' = Q(n, 6) 

I" - Q(n, 7) 

Such correspondence in general program, is obtained by 
introducing some artificial integer variables whose values can 
be fed in the program as lata. Values of these variables are 
different for different set of equations. 

The artificial variable are JJEQ, JJCI, JJCF, NKMP, 

LCH, NIB, NIB, MI, NNV. 

In addition to these variables there are others like 
SST , M, NR, H, EB, JA. 

Some of these variables are double subscripted, right hand 
side subscript showing the set of equations to which these varia- 
bles belong. 

Now a short description of these variables are given b^tlow.. 

JJEQ is a double subscripted variable which is given 
value 0 or 1. If 

JJEQ (n, m) = 1 

it denotes highest derivative of dependent variable occurring 
in one of the equations of m^ set is the n^ dependent variable 
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Hence corresponding to equations (42 ) 

JJEQ (4,2) = 1, JJIQ (7, 2) = 1 
JJEQ (4, 2) denotes occurance of G'" and JJEQ(7, 2) denotes that 
of T" in the second set of equations. 

JJCI is a double subscripted ■variable which tolls the 
program what ■ the missing starting values are.!' ike JT£Q,JJCI 
is given values 0 or 1. Now if 

JJCI (n, m) = 1 

it denotes initial value of the n^* 1 variable is missing' in 

set of equations. Hence, for equation (42), we have 

JJCI (3,2) = 1 

JJCI (6,2) = 1 

These correspond to G" and T' whose initial values are not know 
JJCE is also a double subscripted variable denoting the 
variables which have prescribed values at the end point. If 
JJCE (n, m) = 1 it denotes n™ variable of the m^* 1 set of equa- 
tion, has a prescribed value at the end point. Eor eq. (4 

we have 

JJCE (2,2) - 1 
JJCI (5,2) = 1 

These variables take; the value 1 under the conc’-’tions mentions 
above, otherwise they take the value zero. 

LwH is a subscripted variables. .It can take the value 
0, 1, 2, 3. If LCH (m) equal to zero, m^ set of equation will 
be calculated as will be explained while describing subroutine B 
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If ICH (m) equal to 1 , the set will he calculated by PDM. 

If LCH(m) equal to 2, the set will be solved directly after 

obtaining approximate starting values by PDM. If LCD (m) equal 
to 3, set will be solved directly and values obtained in sol* 

ing the previous set of equations will also be used in solving 
this (m^b) set. 

NIB (m) and NBB (m) gives respectively the number of 
initial and final boundary condition given for the set of 

equation. Here, their values are 3 and 2 respectively for all 
sets of equations. 

NKMP is the total number of sets of equations to be 
solved. Hence here NKMP is equal to 3 ( as three sets of equatio: 
are to be solved here) . 

With these artificial variables, program, reads the identf 
of the equation. Other variables like H,NFV,NR,SST ,EEN supply 
the actual data like range of independent variable t values given 
at the boundary points, etc. 

A short description of these variables are given below 

,NN¥(m) denotes the number of dependent variable as will 

be interpretated by the program in i the m^ set of equation, i.e. 

it is equal to the sum of the maximum order of derivatives of 

different dependent variable occurring in the m*'* 1 set of equa- 

"til 

tions and number of actual dependent variable in the m set of 
equation. 
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Here for the second set of equations 

maximum order of the derivative of G = 3 
maximum order of the derivative of T r 2 
no. of dependent variahles{G and O' 1 ) is 2. 

Hence 

MV (2) =3 + 2 + 2 = 7. 

HR = No. of steps into which domain of independent 
variable is divided. 

H = step size of the independent variable. 

SSI and I'M are double subscripted variables corresponding to the 
boundary conditions given at initial and end points. Here G(0) s 
G'(0), T(0) are given for the second set of equations we assigns 
these values to SST(l, 2) SST(2 5 2), SST(3, 2) respectively. 
Similarly the end point G" (°°) and T (°°) are given, hence we 
assign these values to FM(1, 2) ahd FM (2, 2) respectively. 

First initial value corresponding to the order in which 
different dependent ( variables are arranged (as interpreted by 
program) is assigned to SST (_ 2 . , 2 ), the second to SST (2, 2) 
and so on. The order in which the dependent variable, is arranged 

G, G* , G", G m . T, T 1 , T" 

Now first initial condition corresponds to G, second to G f , third 

SST( 1 , 2) = G(0) = 0 
S ST ( 2 , 2) = G'(0)= 0 
S3T(3, 2) = T' (0)= 0 


to T. Hence 
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Similarly, 

m(i) = a" (oo) = o 

FM(2) = T (oo) =0 

Vrrir’ le PI controls the accuracy to which boundary conditions to 
be satisfied. If MI is given a vain -? 3, min: mum up to third place 
of decimal, boundary condition will be satisfied. 

5 . 3 Subroutines 

There are eight subroutines in this program. They are as 
given below. 

MIST -This subroutine along with subroutines RKM solves the 
equation by PDM. This is effectei .jor the second set of equatio 
by putting LCH(2) = 1 

ANT - This subroutine using the values of JJEQ, JJCI, JJCP, etc. 
defines some other variables like JQ, JO, Ietc. Values assigned 
to these variables are used in the main program and subroutine 
RKM. 

HIX - This subroutine along with RKM, solves the original equa- 
tion directly. Here original equation is the third set. Hence, 
we equate LCH (3) to two. Hue to this value assigned to 1CH 
for the third set this set will be solved directly utilizing 
starting initial condition as obtained solving second set by PPM. 
RKM - This is the most important subroutine which solves a parti- 
cular set of equations depending on the value attained by KMP in 
the main program, KffiP denotes the number of the' set ito be solved. 
When in the main program KMP attains the value 2, 
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second set of equation is solved by the program. As mentioned 
earlier, this subroutine iterates the result to the desired 
accuracy controlled by the value assigned to MI . 

ALEQ - This is a subroutine which solves linear algebraic equatior 
This subroutine is called from BKM for getting new interpolated 
initial values.. 

TICK - This is a subroutine which calculate derivatives one order 
higher. In BKM, this is called to obtain higher derivative from 
lower ones. 

PUNS ON - This is the subroutine where different sets of equation 
are fed in the following way. Let us consider second set of equa- 
tions which can be written as 

q <ii = _ ( ff ’ + (s-f )|3+B(gf ,L HfG"+(T~2f ' G T ) j3)-f n +G"( 1-B) ) 

T” = - ( fs f Pr +B G s’ ?r + BfT’ Pr + T ' - BT'-s') 

Now this can be written in terms of subscripted variables as 

FOMULA = - (F(1)*P(2)+(F(5)-F(2)*F(2))*P+B*(Q(1)*P(3)+P(1)*€ 

+ (Q(5) - 2*F(2)*Q(2)) * p - F(3)+Q(3)(1-B) ) 
where G' M is replaced by a variable FOMJIA 

variable F corresponds to f and s and Q corresponds to G and a., 
subscript of F and Q are decided bjr the same rule as mentioned in 
case of subscript m of double subscripted variable F(n, m). While 
writing down equations in FUNSON subroutine we should remember 
variable like G, T, G' etc. which are actually being evaluated 
in BKM by Bunge-Kutta Method is denoted by Q and other variables 
like f, f 1 , s etc. which has already been calculated and is being 
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utilised for calculating G-, T, G-’ etc. or Q are denoted "by S'* 

The statement hearing formula given above is assigned a statement 
number given by 

I = KMP*100+M 

where M is equal to the position of a equation in the KMP^^ 1 set 
of equations. The above equation is the first equation of second 
set. Hence, 

N = 2*100+1 = 201 

Hence, equation (44) is transferred to HUTSON subroutine as 
( for (3 = 1). 

201 T OMUL A = _ ( p( -| )*f( 2)+F( 5)-F( 2)*F( 2) )+B*(G( 1 ) *F(3 )+F( 1 ) *Q(3) 

+ (Q(5)-2*F(2)*Q(2) )*1-P(3)+Q(3)*(1-B)) 

Similarly, equation (45) is transferred as (for Pr = 1 ) 

202 F OMUL A. = - (P(1)*F(6)+B*Q(1)*F(6)+B*F(1)*Q(6)+Q(6) 

-B*Q(6 ) - F(6 ) ) 

When the third set \s to be’ solved, corresponding equations 
are transferred ' to the FUNS Oil- . subroutine as (for (3 = 1, 

Pr = 1) 

501 JOMULA = - (Q(1 )*Q(3)+Q(2)**2)) 

502 FOMULA = - Q(1)*Q<6) 

BST - If solution for the initial value of the parameter can be 
obtained directly as in the case of equation (41) where solution 
is obtainable in a closed form as given by (41a), then the formulc 
giving the solution for the initial value of parameter is fed in 
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the subroutine BST. This is echioved by giving the value equal to 
zero to the variable LCH for the set of the equation for which the 
solution is obtainable in a closed form. Here equation (41) is the 
first set. Hence we put " 

LCH (1) = 0. 

5.4 Some comments on the General Program? 

To use the general program, most economically and effective!' 
it is suggested that first the problem, should be solved by PDM 
taking large ’step size of the parameter. After obtaining the 
solution at the desired value of the parameter, we can solve the 
original non-linear equations directly with starting values obtaine 
by PDM in the first step. As the assumed initial values are very- 
near the exact valu° , solution will converge to the true solution 
of the desired accuracy in a few iterations. Solution obtained by 
PDM by taking a step size of p equal to .25 and the solution obtain* 
by solving the parent equations directly utilizing the initial boun* 
ary values as obtained by PDM, are given in a tabular form on the 
next page. 

In the main preplan, double subscript variable AD always 
denotes the variable whose values are to evaluated like (f, f ’ , s,s 
etc. The other double subscripted variable is Q which is the vari- 
able whose values are known. In case of PDM., Q denotes variables 
like (df/cip), (df'/dB) etc. 
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Table 1(a) 


ET Li'.}) 

N* J / 

f 

s 

O 

• 

o 

0.000 

1 . 5000 

0. 5 

0.1552 

1.352 

1.0 

0.5171 

1.217 

1.5 

0.9722 

1.113 

ro 

* 

o 

1.460 

1.048 

2.5 

1 . 956 

1.017 

3 . 0 

2.453 

1.003 

3.5 

2.950 

1.001 

4.0 

3.448 

1.001 

4.5 

3.946 

1.001 

5.0 

4.445 

1.001 

5.5 

4.944 

1.000 

6.0 

5.443 

1.000 

6.5 

5.942 

1.000 

7.0 

6 . 442 

1.000 

7.5 

6.942 

1.000 


T Solution obtained by solving 
equation (42) by PDM taking 
B = 0.25 for (3=0,?r=1 


Table 1(b) 


ET A (y'| ) 

f 

s 

0.0 

0.000 

1.500 

0.5 

0. 1565 

1.352 

1.0 

0.5213 

1.217 

1.5 

0.9800 

. 1.112 

2.0 

1.472 

1.048 

2.5 

1.971 

1.017 

3.0 

2.472 

1.005 

3.5 

2.972 

1.001 

4.0 

3.472 

1.000 

4.5 

3.972 

1.000 

5.0 

4.472 

1.000 

5.5 

4.972 

1,000 

6 . 0 

5 . 472 

1.000 

6.5 ' 

5.972 

1.000 

7.0 

6 . 472 

1.000 

7.5 

6.972 

1 . 000 


Solution obtained by directly 
solving the parent nonlinear 
equations staffing value of whic> 
is obtained by solving equation 
(42) by PDM. Here also p=0,Pr=1 . 
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The program can also he utilised to solve more than one 
set of equations where results of one set of equations is to he 
used in tKe;. solution of the next set as is required in solving 
second order boundary layer equations where values from the first 
order is required to solve the second order equations. In this 
case Q can he given values corresponding to the first order 
solutions. This is achieved in the program by giving the value 
3 to the variable LCH for the set of equation corresponding to 
second order boundary equations. 

General program is given in the Appendix 3. This program 
can solve a set of equations containing 25 (computer) dependent 
variables. On IBM 7044 computer,, some 4300 memory cores remain 
unused. Hence by changing dimension,, this program, can handle 
equations with 35 dependent variable. Maximum number of boundary 
condition that can be satisfied at the end point is 10. Time take 
by this program on IBM 7044 to solve equation (42) by PDM and then 
to solve equation (43) directly is 6 minutes 53 seconds. 



DISCUSSION 


Results obtained for Prandtl number equal to unity tally 
quite fairly with We given by Cohen and Reshotko (Ref, 4). 
General trend, of the result ip the same as siven in f\\*»t reference 
(4). Maximum deviation of our result from that given in that 
reference is ,5 percent. This difference can still he reduced 
hy t akin a shorter step size of parameter |3 and that of the 
independent variable , Prom, the result obtained for Prandtl 
numher eual to unitv, we can hope results obtained for other 
Prandtl number are also- reliable so far as the solution of 
equation (15) and ( 16 ) is concerned. 

For Prandtl number not equal to unity, solution obtained, 
in general is the complementary function of equation (15) hut 
gives the near exact solution for low mach number flows. For 
high mach number (Iff - °°) right hand side of equation ( 1 3 ) 
approaches a constant value equal to two. , With slight modifi- 
cation in programme results can he obtained for high mach number. 
Thus solution obtained for Prandtl number not equal to one may 
not give exact quantitative value in all cases. Nonetheless, 
values were calculated for Prandtl number not equal to unity 
as it gives a qualitative understanding of boundary layer. 

Our main aim was to solve the compressible boundary 
layer equation by parametric differentiation method and we shew 
that this method can be effectively applied to solve the problem 
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APPENDIX -l 



"TAgte / 

SIFT A = 0.0000 PRANDTLE NO. = 1.0000 

vti.cniry gradient at wall = 0.4696 

T » M p r R A T 1 J R F GRADIENT AT WALL = 0.4695 


DISPLACEMENT THICKNESS * 

-0. 

0001 


MOMENTUM 

THICKNE 

ss = 0 

.4696 



THE R M A 1. 

THICKNE 5 

S = -1. 

2169 



REYNflL D 

ANALOGY 

PARAMETER = 

2.0004 


ETA 


F 


DF 


0. 


0. 


o. i42.ee- 

07 


(. j 

0.5364!=- 

CL 

G .23425= 

00 

■ . ] C ; ; < - 

T i 

0 .233 0 c 

00 

0.4606E 

00 

. l v: of 

1 

C.5150E 

00 

0.66I5E 

UO 

O.REOn: 

0! 

0.886 BE 

00 

0.81 6 7 h 

00 

0. 0 5 

■:i 

0.1322!- 

01 

0.91 68 E 

00 

. .-.i o o 1 ’ r 

: li 

0 . 1 79 6 E 

Oi 

0.9 6 9 IE 

00 

0.3‘>L9 r 

- Oi 

0 . 2 2 8 6 E 

01 

0 .9907 E 

00 

1 1 » 4 ('; i ) r :• ! 

f \ 

0.27840 

01 

0.9978 £ 

00 

0 . 4 ‘ ‘ 0 ' 

■ 01 

G.32R3E 

0 1 

0.9996E 

00 

o . T C -j o r 

0 1 

0.37B3E 

01 

0.9 9 99 E 

00 

j, 6 50 05 

i. 

0.42 8 3 E 

01 

0 . 1000F. 

01 

0. CO GO! 

: oi 

0.478 3 E 

Oi 

0. 1000 E 

01 

0 .6 ':.0l ! 

'■ u. 

0.5283E 

01 

O.1O0OE 

01 

0 . 7 1 o n i 

1: Cl 

0.S783E 

01 

0.1000E 

01 

0,7 50 0- 

t Cl 

0.6283E 

01 

0.1000E 

01 


s w 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0.0000 


s 

. 2 342 E 00 
.4606E 00 
.66 14F 00 
. 8 1 6 7 E 00 
. 9 1 6 B E 00 
. 9690E 00 
.9 90*7 E 0 0 
. 9 9 7 8 E 00 
. 9996E 00 
.9999E 00 
i.LGOOE 01 
t. 1000E 01 
I.IOOOE 01 
)* 1000E Cl 
I.IOOOE 01 


-TAgLE 2- 


BFTA = l . 0000 PRANDTLE NO. = 1.0000 SW = 0.0000 

VELOCITY GRADIENT AT V. ALL = 0.648,6 

TEMPERATURE GRADIENT AT WALL = 0.5065 

i) I S PLACEMENT THICKNESS = -0.1575 

MOMENTUM THICKNESS = 0.4036 

THERMAL THICKNESS = -1.1374 


6YNCLD 

ANAL DGY 

PARAMETER = 

2.5612 




CIA 


F 


OF 


s 


0 . 


0. 


0 * ' f> ’) -.1 *“ 

07 

0. 


V. , s; T ( ; , 0 

i 

0 . 7987 fc- 

01 

0 .3149 E 

CO 

0.2524E 

00 

*' * 1 C • ; r • i 

i 

C . 5 0 6 7 E 

GO 

0.6817E 

CO 

0.4936E 

00 

iu 1 r vA-\'. 

h i 

G . 6 4 9 ^ E 

00 

0. 7 7 6 2 E 

00 

0.6998E 

00 

■■ . * / '■ - ' / 1 

01 

0.1071E 

Cl 

0.3 97 IE 

00 

0.8492E 

00 

r..?*vu > 

01 

0.1537E 

0 1 

C. 9602 6 

00 

0.9374E 

00 

•y t r \ v 

■# 

0 1 

C.2025E 

01 

0.9873E 

00 

0.9789E 

00 

0, r.i 

01 

0.252 IE 

01 

C.9967E 

00 

0.9943E 

00 

. . */» r: • 

0 1 

0.2020E 

01 

0 * 9 9 9 3 E 

00 

0.998 EE 

DO 

o • / ; y.;Or 

Cl 

0 . 3 5 2 C £ 

01 

0.9999 E 

CO 

0.9998E 

00 

<• *0 t ; 

.* # ■* 5 • ■? 

01 

0.4 02 OF 

01 

0. 1000E 

01 

0. ICOOE 

01 

r. % *■; !'■,'}(, i 

0 i 

0.4 5P0E 

01 

0. I 0006 

01 

0. 1000E 

01 

l . 600 or 

Cl 

0.5020E 

01 

0. 1000E 

01 

0. 1000 E 

01 

r, , 6 yOGi„ 

C 1 

0.552QE 

01 

o. icooe 

01 

C.tOOOE 

01 

0.7 CO OF 

■ Cl 

0 * 6020E 

01 

0. ICOOE 

01 

0. ICOOE 

01 

0 . 7 6 0 0 E 

. Cl 

0.6520E 

01 

C, 10006 

01 

0.1000E 

01 
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TA&Le 3 

P'FTA = 3.0000 PRANDTLE MO. = 1.0000 SW = 0.0000 

\/lLt:01TY GRADIENT AT WALL = 0.7384 

T fc M PEft ATUR6 GRADIENT AT WALL = 0.5205 


:> i J L AC F 

MH u thickness = 

-0. 

2060 




M tz fs.J J I J,M 

THICKN6 

SS = 0 

.3839 





r R M A L 

THICKNESS = -1. 

1109 





ynhl n 

analogy 

PARAMETER = 

2.8374 




E T A 


F 


DF 


s 


# 


0. 


-0.3332E- 

•07 

0. 


r" ^ r* /*\ 

c o 

0.0001”- 

„ . 

0.3517E 

00 

0.2693E 

00 

•% :' r • , ’ 1 . * 
* i ' . ' .* ! ' " 

i: i 

0. 3385E 

r , 

s; 

0.62706 

00 

0.5 0 5 6 E 

00 

, 1, G\: *.r. 

0 1, 

0 . 7 0 1 4 Ei 

0 0 

0.30996 

00 

0.7131E 

00 

r ? ' , • ■ 

r ] 

0. ] 1361= 

01 

0.9 1566 

00 

0.8596E 

GO 

.#? * 

r. ] 

0. 1 60 8 E 

01 

0.9682E 

00 

0 .94346 

00 

. :• 0 ■’ : ' r 

G 1. 

0.2099E 

01 

0 .99006 

00 

0.9F15E 

00 

. 7 

C 1 

0 . 2 5 9 6 £ 

0 1 

C. 99750 

00 

0.9952E 

00 

. A*. cor 

C .!. ■ 

C. 30956 

0 i 

0.99956 

00 

0.9990E 

00 

. 4 r - f; ir 

\ 

0.35956 

01 

0.9999E 

00 

0.9998E 

DO 

♦ "1 ,* i. " 

G ], 

0.40956 . 

01 

0.1000E 

01 

0 . 1000 E 

01 

■■is O, , ' ! J I 

01 

0 .45956 

01 

0.10006 

01 

0. LOO 06 

01 

. oooo c 

0 1 

0.5095E 

01 

O.IOOOE 

01 

O.iOOOE 

01 

•. 6*500*: 

i; 1 

0.5595E 

01 

0. 1000E 

01 

0. 1000E 

01 

. /oo or 

c 1. 

0.60936 

01 

0.1000E 

01 

O.IOOOE 

01 

. 7 5 ; i 0 E 

n. 

0.65 9 5 £ 

01 

0.10006 

01 

0. 10006 

01 



TABLE 4 

CCTA = 0.0000 PRANDTLF. NO. = 1.0000 SW = 0.8000 

VELOCITY GRADIENT AT WALL = 0.4696 

TEMPERATURE GRADIENT AT WALL = 0.0939 


DISPLACEMENT THICKNESS = 0.9734 

MOMENTUM THICKNESS = 0.4696 


T HER MAI 

THI. CKisiES 

S = -0. 

2434 





RE YNCLD 

ANALOGY 

PARAMETER = 

2.0004 




ETA 


F 


OF 


s 


0 . 


0. 


0.14286- 

07 

0.8000E 

00 

c . scor.r 

C 0 

C . 5864 E- 

•01 

0. 2 34 26 

DO . 

0.34686 

00 

V . ] CEDE 

• , 1 

0 .233 0 E 

0 0 

0 *46066 

GO 

0.89216 

00 

. 1 SCO!;' 

c 1 

0 . 5 1 5 0 E 

00 

0 . 6 6 1 5 E 

00 

G.9323E 

00 

0.?Cu,.T 

C 1 

0.8868E 

0 0 

. 0.81676 

00 

0.9633E 

00 

U :i.l 

0 1 

0. 13? 2 E 

01 

0.91686 

00 

0.9834E 

00 

’ ' % T • ? J ■ < \ 

Cl 

0, 1, 7 9 6 E 

Oi 

0.96916 

00 

0.99386 

00 

. 3CGU ! ." 

Cl 

0.2286 E 

0 1 

0.990 76 • 

00 

0.998 IE 

00 

C .iCE 

0 1 

0.2784E 

01 

0 .99786 

00 

0 . 9996 E 

00 

f;, . A 90 DV 

0 1 

0.3283E 

01 

0.99966 

00 

G.9999E 

00 

< : . 9 C E- D 1 : 

a 

0.3783E 

Oi 

0.9999E 

00 

0. 1000E 

01 

* / « S r * 0 oc 

0 1 

0.4283E 

01 

0. 1000E 

01 

O.iOOOE 

01 

C . 6 0 0 0 T" 

CL 

0.47B3E 

OI 

O.IOOOE 

01 

0.1000E 

01 

U.6C0CC 

Cl 

0.52836 

Oi 

0.10006 

01 

O.IOOOE 

01 

I/# TOC-D-C 

Cl 

0.5783E 

01 

0. 10006 

01 

O.IOOOE 

01 

C.TSoOC 

01 

Q .62836 

01 

O.IOOOE 

01 

O.IOOOE 

01 



B f T A = 

TAgLS 

I • 0000 PR A NOT LE NO 

s~ 

. = 1.0000 

SW = 0. 

,8000 

vti.nc i ry 

GRADIENT 'AT WALL = 

1.1234 




TF.MPFR AT 

URL GRADIENT AT 

WALL = 

0.1119 




DISPLACE 

ML NT THICKNESS = 

0.4950 




MOMENTUM 

THICKNESS = C 

1 . 3 1 6 3 





THERMAL 

THICKNESS = -0. 

2 08 8 





REYNOLD 

ANALOGY PARAMETER = 

4.0165 




ETA 

F 


OF 


s 


0. 

0. 


-0.9810E- 

07 

0.8000E 

00 

0.5000E 

CO Q.1238E 

c c 

0 *4622E 

00 

0.8556E 

00 

u.i'cvor. 

01 0.4320E 

0 0 

0.7442F 

00 

0.9076E 

00 

W. ! B B Of 

il 0.845 BE 

on 

0.8934E 

00 

0.9494E 

00 

• 2 C l . / 0 >; 

Cl 0.1312E 

01 ' 

0.9616E 

00 

0.9770E 

00 

u . ? V,; or 

Cl 0.18 0 1 E 

01 

0.988 IE 

00 

0.9914E 

00 

0 * '^OUOt 

Cl 0.22 07 F 

01 

0.996 9 C 

00 

0.9974E 

00 

0. 5B r: 

01 0.2796E 

01 

0.9993E 

00 

0.9994E 

00 

0 » f\ t ‘ 0 0 B 

Cl C . 32 96 F 

01 

0.9999E 

00 

0.9999E 

00 

{■ * A '0.>Oh 

Cl 0.3796F 

01 

0. 1000 E 

01 

O.IOOOE 

01 

0 • ’> 0 u 0 ; ’ 

Cl 0.4296E 

Oi 

0. 1000E 

01 

O.IOOOE 

01 

■J * t) ^ •'./ ! 

01 Q.4796E 

01 

'o.ioooo 

01 

O.IOOOE 

01 

0 * C U 0 ■ ' b 
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P> I- T A = 0.0000 PRANDTLE NO. = 0.3000 SW = 1.5000 
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\\ r T A * 0.0000 "p'UNGTti HO. = 9.0000 SW = 1.0000 

Vi'.iCf.I ry GRADIENT AT WALL = 0.4696 

T i: M P LU AT UR!' GRADIENT AT WALL = 0.0000 

n I P L A C L MENT THICKNESS = 1.2168 


H )*'!. N I UM THICKNESS * 0.4696 

THfPJOH. THICKNESS * -0.0000 

i: T A F 

* » 

0 . !5 364E-0 1 
• . ] r • S C.2330E CO 

. ! 0 .51 r > Z F 0 0 

: 1 0 , wp.68f; 00 

- r? 0.1 32 2 E 01 

, , : 0.179&E 01 

> , •;». 0.2266c 01 

. . :i : •' | 0.2 7P4F 01 

1 0.32 8 31; 01 . 

mr < I 0,37835 01 

,,9,:.,=. . i 0.42B5E 01 

;. » e,C :,"A : l . Q.4783E 01 

4^ OOF L-.l Q.5283E 01 

7 ,;- {i v‘ cl ' 0.57835 Qt 
.. , ( \ . 0.6283F 01 ■ 


OF 

s 


0.14285-07 

O.IOOOE 

01 

0.23425 00 

O.IOOOE 

01 

0.46065 00 

0. 1000 E 

01 

0 . 6 6 1 5 E 00 

0. 100 OF 

01 

0.8167E 00 

O.IOOOE 

01 

0.9168E 00 

0. 1G0QE 

01 

0.9 69 IF 00 

O.IOOOE 

01 

0.9907E 00 

O.IOOOE 

01 

0.99 7 8 E 00 

O.IOOOE 

01 

0.9996E 00 

O.IOOOE 

01 

0.9999F 00 

O.IOOOE 

01 

O.IOOOE 01 

O.IOOOE 

01 

0.1OOOE 01 

0*. 1000E 

01 

O.lOGOf- 01 

0. 1000E 

01 

O.IOOOE 01 

O.IOOOE 

01 

O.IOOOE 01 

O.IOOOE 

01 



t ^ $ LB 

or: r “ " 0. 9000 PRANDTLE NO. = 9.0000 SW = 1.5000 

V> UK I IV GRADIENT AT WALL = 0.4696 

TTWRATOHE GRADIENT AT WALL = -0.4965 

niSPi \Gf 6 N T THICKNESS * 1.5020 

NDMFDTtiO THICKNESS = 0.4696 


f M /t ) 

THICKNE 

ss = 0 

.2852 





YMDL 0 

ANALOGY 

p * *"’* r. * t 

7 f « : = 

0.9458 




i. T A 


F 


DF 


s 




0. 


0.14296- 

-07 

0.15006 

03 



9. 5564 E 

-01 

9.23426 

00 

0. 12 57E 

01 

1 ; ’ 

*• 1 

: i 

0.2330!= 

00 

0.46Q66 

00 

0.1075E 

01 

. • < : ■' 

! ■ l 

0.51 506 

0 0 

0.6615c 

00 

0. 1008 E 

01 

♦ ; 1 • 

I 

0.P868H 

00 

0.81676 

00 

0. 1000E 

01 


1 

0. 1322F 

01 

0.9 1686 

00 

0. 10006 

01 

® ■! .. ' 

i. 1 

0.1 7 76 5 

01 

0.969 IE 

00 

0. 1000E 

01 

« / . ; 

. 1 

0.22 861: 

01 

0.99076 

00 

0. 1000E 

01 

' ■■ 

, 1 

0.2 7846 

0 i 

Q.9978E 

00 

0.10006 

01 

* ’i ’ f 1 : ’ 

: . 1 

9. 3 2836 

01 

0. 9996 E 

GO 

UJ 

o 

o 

o 

* * 

o 

01 

# " % ’■ . . ' : 

• 1 

0. 37832 

01 

0.79996 

00 

0.10006 

01 

«F i.*< ■ ' ’4' . 

U l 

0.4 28 IE 

01 

0.10G0E 

01 

. 0. 100 OE 

01 

. 6 : - ■ : ■: f; 


0.4 7510 

01 

0.1 0006 

.01 

o.ioooe 

01 

, (> m i) 

* 

0.6 2 8 3 E 

01 

o.ioooc 

01 

0.10006 

01 


l 

0.87836 

01 

O.'iOOOE 

01 

0.10006 

01 

. rvUO'’ 

.•a. 

0.6283E 

01 

0.10006 

01 

■■ Qm 1QQ0E 

01 



£7 
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H r T A - 

TABLE 1ST 
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4;; 01 0.7044E 01 
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i' r T A * “f!.0 625 RRANDTLE. NO* = 1.0000 SW = 0.8000 

Vi ICC 1 TY OR All ] ENT AT WALL ■= 0.3952 

Tl-MPf 4 i! •- GRADIENT AT WALL = 0.0908 

0 I SIM A C I ’O f iT n 3 1 .0718 


m Ml N ! I, THICKNESS = 0.4929 

T \ i ! 9 MAI THICKNESS = -0.2507 
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ANALOGY 

>■» A P ARFT'' 
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A At. 

0.517 75 

01 
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A 0 1 
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TA 62 -/' zj 

ft F T A -0.1250 ^RANOT'LC ' WO. *. 1.0000 SW = 0.8000 

vrtnOilY CHAO I «r NT AT WALL a 0.3049 
T f; f" i 1 1 ' , ’ A H i . ’ UR At) IFNT AT WALL = 0.0865 

rm>n.ACJ ms : ?u r,i TCK^rtis = 1.2141 

MUMr jm: THICKNESS * 0.5217 

THFWMAt THICKNESS * -0.2617 


YNC1 f * 

1 ANALOGY 

0 A, p- ft !/: F r £ R = 

f 

1.4102 

OF 

s 

* 


— 0 * 

("> • 4 ' ) 1 9 F — 0 1 

-0.1019E-07 

0.16406 00 

O.bCOOS 00 

0.8432E 00 

!: | '» 

, ■ 1 

0.1 6 048 00 

0.3506E 00 

0 . 6 8 5 3 E CO 

* i i 

• ' : !. 

0.3 9? Oil 0 0 

0 • 6 4 3 1 E 0 0 
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• j ; • 

• . . . | 
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■ tt / " ‘ ; 
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0.97776 00 
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-.4! , 

' '1 ’i 

ij. 25 26b 01 
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; ' A . ■ 

. .* r ( t | 
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i ; * It c 

1 i ’ V Cl 
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U * ** - 
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. 0 * 4024 E, .0 1 
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0. 10006 01 
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5 ,: v A! 
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‘ 0.1000E 01 
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cl. 
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0.1 0006 01 

O.IOOOE 01 

Q * 7 Vi 

:voc cl 

, ■ . 0 • 5 5 2 4 E CM . 

0 .10006' 01 

O.IOOOE 01 

* ,7V 

Dnr ■ c-l ... 

. .0., 60246 01 

0* LQOOE 01 

O.IOOOE 01 



■ 2 .^ 

= 0.3000 SW = 0.8000 

0.4149 
= 0.0597 


0.9015 


~T AgL e 

BETA = -0.0500 PRANDTLE NO 

VELOCITY GRADIENT AT WALL = 

T t-. MPERATURE GRADIENT AT WALL 
H I $ PLACEMENT THICK 
MOMENTUM THICKNESS = 0.4840 

THERMAL THICKNESS = -0.3866 

REYNOLD ANALOGY PARAMETER = 
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VLIUOiTY GRADIENT AT WALL = 0.3814 

rtNPCl'.A TORE GRADIENT AT WALL = -0.0000 
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MOMENTUM THICKNESS = 0.4962 

T HER MAI THICKNESS = -0.0000 
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p E T A = -0.0625 PRANOTLE NO-. = 9.0000 SW = 0.8000 

\l F I. CC I T Y OR AH IF NT AT WALL = 0.3888 

TEMPERATURE GRADIENT AT WALL = 0.1889 

D 1 S PLACEMENT THICKNESS = 1.2151 


MOMENTUM THICKNESS = 0.4956 
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FORTRAN SOURCE LIST 

’ sn snu.Rcr statement . ■ 

0 T. I8FTC ,'•• ■ 

1 DH-cNSIlJM F{ ? 1 0 ) ,F1> (210) , F? (210) ,F3(210) , S(?10) , SU210) ,S2( 21 

3 S 3 ( .? 1 0 } * G 1 ( ? 1 0 ) » G 2 ( ? 1 0 ) » G 3 ( 2 1 0 ) » T 1 [ 2 1 0 ) * T 2 ( 2 1 0 ) » T 3 ( 2 1 0 ) » 

?■ A h T A ( ? 1 0 ) » D F ( ? 1 0 ) , B 0 (210) 

3 " ( o) ! b) t A A ( t> } , L IE (210) , AFPI2101 , SWW { 6 1 

2 LOGICAL LAO 

3 4 Ei 'Kka T ( /?2X » 2 7HVELCC I T Y GRADIENT AT WALL = t F9.4) 

4 3 FCJRp ATI/ 2 r X » 3 0 H T F M ° T 3 A T U R E GRADIENT AT WALL =,F9.4) 

10 FOR fcA T ( /? 5 X , 3HE T A , 1 2 X , 1 HF , 14 X , 2HD F , 13X,1HS/2X) 

6 f: F C R ? : I { / 3 0 X , 2 4 h D I S P L A C F M IF' ! T THICKNESS =,F9.4) 

7 9 r-OR' ; AT(HU ,2CX,6Hr;FTA =, F9.4 ,3X t 1 4H -‘RANDTLE NO. =,F9.4 ,3X, 

] ~ » r • a ) 

I r 7 Fl'l’i A T { / 2 0 X T 2 0 H M 0 7 E N T ! v. THICKNESS = , F9.4) 

I I 3 FORMA r ( /?CX , 19HTHCRMAI. THICKNESS =,F9.4) 

12 99 Fi 1 ft •' i a T ( / 2 OX , 2 7HP. E V MC L L ANALOGY PARAMETER =,F9.4> 

13 II FOR n A T ( / .? o X»F11.4»4X»H11.4»4XfFll.4»4X»61l.A) 

1.4 a EAR 9 99 ,H,DE 

1.3 'A 14 FORMA T { 4E 15 • 5 } 

14 999 FuRHAT(F4.2 t F4.2) ■ 

17 777 F0kMA7(9E1 2.4) 

20 55 b r p 3 Y AT { F 3 . ? » 1 0 X » E la . 5, 10X » El 5 . 5 » 10X , I b ) 

21 419 FORMAT (3215. 5) 

?.?. CALL FLU.N ( 2 SCO 5 

23 KST=l 

24 S W W ( 1 ) = 0 • 2 

2*5 SWW 1-2 1*0.6 

26 " SWW 131*1.5 

27 SWW(4)=I„ 

30 SwW(5)*2. 

31 IKK *4 ■ 

33 PR =5 . 

3 3' 147 HO 3 3 1 M$T= 3 1 IKK . 

34 J A = 2 1 

35 ' NA* 151 

36 4=1. 

27 ■ AETAU)'»G. 

40 . 00115.1=1 » 2 ... : : 

■41/ . J= I 

42 v., . S'= 0 . 

4 3 7=0. 

■. 4 4 ' W~ A 

45 ETA=0. 

46 N = 1 

4 7 F ( N ) ~0 . ' 

50 222 CKO=H*Z : , ■ ' 

51 OL0=H*W 

52 OM0=-H*<F(N)4W) . • /' 

53 0K1=H*( Z+DLO/2. ) 

54 DL 1 =H’T ( W+UMO/2* ) ■ 

55 cf/ 1 =-h* C C *" C N 1 +CK0/2 . 1 * I W+0M0/.2 . 1 1 

5 5 0 K. .7 = H ’!* ( Z + L 1 / 2 . ) . ' v .v. ,• ' ■ •/; ; • 

0 7 OL? = h* ( h+DMI/2. ) , 

t.O T;M2 = -H* ! { f ( H) +0K172. ) *( W + 0M1/2 . ) 1 

1 ] i ) k. 3 = H 7 { / + 0 !. 2 5 

f,P P L 3 = H * ( V'i 4 D 0 2 ) 


1 /, 


fortran source list 

•? r '•* 1 

1 *5 < 


sour , cl s T o t e f ' m t 

6 7 


r.;3«-H* ( ( F ( 0) + PK2 ) * ( W+DM2 ) ) 

(A 


F ! < v + 1 ) =F ( M) +1 . /6. * { DK ;■•+? ,*DKl+2. *DK.2 + DK3) 

6 6 


l = 7 + 1 . / 6 . * ( D 10+ 2. * D L 1 + ? . * 0 L2 + DL3 ) 

6 6 


w = v + 1 . / 6 . * ( 0 f/ 0 + 2. * 0 r 1 ± ? . * ly - 2 + H v ) 

hi 


IT I M- " ) 11 3 ? 1 1 4 » 1 14 

7 f ■ 

1 1 3 

r TA=r TA + H 

'? ] 


‘Vr. 14 4*1 

7 ? 


CL' TV: 222. 

7 3 

1 14 

*■* ( l./Z3**1.5 

7 4 

1 1 6 

r T‘. . T T Ml 1 r 

0.- ^ 1 }. 1 J, ’ *s< l,> 

76 


CALL T ICMHt’NA, F, OIF- ) 

7 / 


J 1 = 1 

1 c o 


5 K = 1 

j J 


MK = 1 

1 7 ? 


S6=5LH(0ST) 

1 i: i 


0=3 .. 

l C* 


SI 1 ) =sw 

J (; r . 

*r i 

r ) f-y 

o s w 

166 


L = 1 

3 t.7 


t:N(l}=o. 

! 1 0 


ro ( i ) =n. 

1 1 i 

K C 

A = 0 

‘ 12 


* * [j c II I ^ 4 f 

1 t 3 


CALI FLUN(IOO) 

i .1.6 


ir- 1 1 . g r . 2 ) a = i . / ? . 

LIT 


IF I 1 ,GT.4 3 A = l. 

1 2 2 


■FLLL*FfL)+01F(LJ*A*H 

1. 2 


DM m = I P+DM U - 1 ) * A ) *H 

1 2 6 


ON ( ! ) s-H4pfI#FLLL* t P + DN ( 1-1 ) *A ) 

125 

2 1 1 

CONTINUE ■ 

1 2 7 


on M* ( 1 . /& . 3 *( DM ( 23 + 2 . *DM (33+2. *DM { 4 3 + 0 M (533 

1 3C 


ODM = ( l«/6. )*( DM ( 2 ) +2 . *l.)N 13 3+2. *DN (43 + DM (5)3 

1.31 


S l L + 1 1 = S ( L ) + DDH 

1 32 


paR+DCN' 

133 


ss-sa+i.) 

1 36 


IF ll.-NA36.6-f 77 » 77 

13? 

46 , '' 

l = L + l 

134-. 


GO T O' 5 F> ■ . .. ■ ■ . 

1 37 


I r ( A is S' ( s { NA ) - 1 . 3- . 00001 ) 92 , 92, 88 

■1.4-Oi ■ 

j’Sv.v 

AAt I I ) = $ ( N A 3 

161 


I F ( I I -2 3 6.9 1 9 C ♦ 90 

i42 ; 

: 'a '9 c 

11=2 

143' 


W = W+0.'5C00 

1 4 4 


NK*fcJK+l 

! 6 5 


CO TO 5.6’ : 


9 9 

IF(F.K.FQ.2)GCT092; 

1 '3 1 

91 

I F { &P S { AA ( 2.) -A A { 1 ) } . LT ..013 GOT 089 

1 56 


W= { A ?. ( 2 3 “ W* A A ( 1 3 +W-i. 3 / ( 4 A ( 2 3-AAl 1 ) 3 

I 3 5 


K K. = K K •+ 1 ... 

1 56 


GOTO? 6 . ' ; . : . 

1 •> 7 


Or; 4 4 4 K= 1 » NA ' .4 

i 60 

4 4 /* 

AH T A ( K+ 1 3 = AfcT.M K 3 +H . 4 


& 


HE GO 04 


•URGE ST". 


FORTRAN SOURCE L I ST 


I F { KS T. IT . 3 ) 00102 

)? DO 10dM-ltNA 
F 1 (M) =F CM) 

SI ( M ) = S ( M ) 

)0 CONTINUE' ■ 

C Al LKUT TA { NA , !>R , FI , s 1 , B , H, G1 , T1 , AZ , BZ , LM ) 

001 01 M=J. ,MA 

F2 ( M ) =F 1 { M) +G1 (M)*DB/2. 

)1 S2 ( M 3 =S 1 ( M ) + T1 ( M ) *DB/2. 

C=B+DB/2. 

C A LL K UT T A C NA * PR , F2 » S2 » C t H t G2 * 12 * AZ » BZ » LM 3 


>Sd 

p ■> 

( 0 3 = 

; r?{pj+G 

2 (M)*DC/ 2 . 


.. 0 6 
■ i d 

). * • ,$ 

1 3 ( M ) 

: ' 1 + ( ' '■ 1 

=S2(M)+ 

T 2 ( M ) £ 0 & / 2 . 


• 1 1 

C A 

1 1 K 

PT T A { NA 

» P R j F 3 } s ? C 

» H » G 2 , T 3 » A Z » 8 

1 1 2 

; i ! 

1 ,V 

L , , * 




l 3 

t C 

M) = 0 

• / 6 . * ( C 

1 ( M 3 +4 . * 02 ( 

t- )+G3(M3 3 + F ( M 

‘ 1 *+ 

i .; ? $ { 

* V - > = <' 

■3/6.* ( T i ( M ) +4 . * T2 { 

M 3 + T 3 ( M } 3 + S ( M 

* I (■: 

IF 

( j-j 

A ) 1 0/ , l 

05,103 


f I 7 

j r t 4 p - 

P 4* P P 




: 2 ■ 3 

*i *** 

J+l 




' ? 1 

i r 

((Kb 

T.Fo. 1 ) 

. AND. ( J.EQ. 

113 3 GOTO? 

2 A 

IF 

( 1 F S 

T.Fn.n 

.AMO. (J.EQ. 

213 3GOTG2 

■ 2 7 

L 

■-".= ! 

I.F0.23 

.OR. ( J.F 0 . 6 ) .OR. (J.EQ. 8 ) 

■ 

1 1 

( ( K S 

I * t: 0 m2) 

.AND. LAG 3GCT0? 

■ "'4 J 

IF 

( ( K $ 

T . Fu, 2 ) 

.AND. ( J.GT. 

8 ) 3G0T0342 

■ .v« 

I F 

( ( K S 

T.F0.33 

.AMD. (J.EQ. 

11 ) 3G0T02 

■4] 

If 

( C K 5 

F.F0.3J 

•AMD. ( J . 20 . 

21 3 3GOTC2 


If 

{ ( K 5 

F.B/.4 3 

.AND. ( J . EG. 

11) 3GOT02 

A 7 

rr : 

t (KS 

T.Fg.4) 

.AND. (J.GT. 

1 1 3 3GOT0342 

■ H ? 

rr 

S {.KS. 

T.F0.5) 

.AND. ( J.FQ. 

11 ) ) GO T 02 

P !:» 

/ A 

IF 

•/** i • 

( ( K S 

T f j 1 

r.so.53 

7 

* A N 0 • ( J ® 0 T * 

1 1 3 3 GO TO 3 42 

4 i \ f 

■6 1 

\ti ■ 

C A 

1 •.) t l., 

L L T | 

C K ( H » MA 

, F , F 1 3 



no 

30C N 

T*1 ,NA 



6 H 

3 ■" : '"3 A 

FP (0 

T 3 - -F 

1 ( N T 3 + S { N T ) 


6 f 

C.A 

I..L S 

Ul-( AFP » ’ 

NA, H.» ZA 3 


t 6 

on 

V'j t. h 

T.*l »MA- 




AFP (NT ) = F1 ( NT 3 *(1. -F I (NT ) 3 

C ALL. SUB ( A PP , NA , H, *6 ) 
CALl$UB(S,NA,H,/.P) 

ZP--/.P-7.5 ' 

NH = 1 * / ( 12 » *H 3 

ZH= ( -25 . * F 3 ( 13+48. «F1 (23-36 
ZR* ( - 20 **,$ ( 1 3+48,*$ (23-36 


, MA , H , P 3 

.*H) ^ " ", 

Pl'( 13+48. *F1 (2 3-36. *F1 (33 + 16. *F 1 ( 4 >-3.*Fl < 5) )*HH 
S (13+48, *$ (23-36.*$ (33+16. *S (43-3, *5 ( 5 


77 

•5 02 

• •' If 
If 

{ HST.Nf.4) Z'Q; 
(M < >t.^NE..4)ZM : 

)0‘> ■ 

■ mi 

3 83 

306 ' 

PK 

mry-fHi pn,sw 

3 07 

PH 

1NT4 ,.ZB: \ 

1 0 

'. PH 

WT5:tM SO:' . 

3.11, 

■■ .■■" '■ fif 

TNT6,Z A 

31? 

"■ fly 

I « TZtZD 

213 ' '■ 


I-!T,8,/P 


WGOCsA 


I su 


S 0 U ft 0 ft 3 T T ?/ ^ jM J 

:* 1 4 


! F ( FST.LT. A ) PH I NT 99 

317 


PRINT IQ 

3 ? 0 


F T A = 0 * 

3 2 1 


LHi 30AN= 1 » NA » lO 

12 ? 


PR I.MT 11 , E TA , F { N) , FI 

3 2 3 

3 0 

A f T'i = F T A+ . 5CCC00 

3 25 

3 ;3 

3 CONTINUE 

3 27 


GO Tru e? 

3 30 

1 0 5 

GO TO 3 31 ■ 

3 3 1 

33 1 

CONTINUE 


•a A ? 

I F ( ! 'ST, F 0.7 ) GCT0321 

> 3 6 


I F ( K S T . t: Li . 3 ) G 0 T Q 3 a i 

ft 3 3 


3r 3 K FI. A) GOTO 3 51 

ft 4 *■■ 


T r ( !< f T ,rr lm q ) r.HTO'jftl 

' 4 7 


OR 

1 3 r ‘ 


i< o T - 7 

1 1 


I R ■- '3 

a n 3 


G ‘ ■ T 0 1 A 7 

> 3 ? 

? ] 

NL=C.l 

i ^ /j 


■< A T = 7 

s b 3 


c g = y . 

*. 3 ( • 


r k k - 

■> s ' 


GFiTrUT 

•5 60 

■* A ? 

* .*» 

jC J ~ ^ ■ 

; 6 1 


r i ■ = . i 

vt 


P j.*, = . 3 

"•t £-i ’■ 


I K * - 3 

■! 6 A 


GOTO! a 7 

■':• o 0 

'1 h ] 

'< S T = A 

■' f;, f; 


1.0 = . 1 

6 1 


'•»«*«. y* - 

7 i ; 


I K K s 3 

-s 7 1 


GO T 01. A 7 

'■?.? 

301 , 

STOP 

, 7 ? 


run 


FORTRAN SOURCE LIST 


'MFGOO'i ' . 

NO MJ-SSAC.CS FOR ABOVE ASSEMBLY 


IBMAP ASSEMBLY 


// 


FORTRAN SOURCE LIST 


STATEi^T 


‘EGCJQ^ 


Si TO- nt:T I '.<2 FL I CK ( H » N A » AF > OIF » 02 F ) 
OP'UINSICN A F ( 2 1 0 ) ,D1F(210) ,D?F(210) 
T SLIP LU: M { SCO) 

\ vA=NA-A 

iVS !\ + ]_ 

H i * = 1 . / ( 

!-!;! = ’ . / ( I 2 . * ( H **2 i ) 

H,j 1 = 1 . / (1 2 .*((?. *H ) **2 ) ) 

FK = ( 1./?. 


OH L-l *4 _ , AC H . 7 v +1 i # a,f ( L + R )-3 .*AF ( L +4 ) ) # HF 

i 4 f 1 » rf* A f { L+ 4 ) ~ 10 * * A F ( L + 5 ) ) 

" Y?=H J 1 * ( 45. *AF {L)-154.*AF( L+2 \ + l] 

1 -IF6.*AF( L + 6) +61 .*AF ( L+8 )-U •-AF ( L *- 5} 

VZf (| } = ( HK* Y2—Y1 ) /HK 1 

Y1-HUMAF(L-?)-8.*AF (L-l I+B.^AFIL+l )~ * .*AF(L+2) ) 

X? =H J ❖ { AF ( L-4 ) — 8 »* AF ( L~ 2 )+8.*AF(L+2 )— i.*AF(L + A 

DlF d ) = !HK«X2-X1)/HK1 F(L) + l6 ^ Ap(L + 1) _ AF(L+2 ) ) 

D?f (t ) = (HK#Y2~Y1 ) /HK1 

DOOta ••<».** *af(L-?|-16.*AF(L-3)+3.*AF(L-4) I 

r vi^[:!'X"r F a!'ui. ! AF««-t.* 2 l*.*AML-2.-i S 6..AML-3» 

l+t,l.*AF(L-4)-10.*AF(L 5) > _ +214 .*AF (L-4)-156.*AF ( L-6) 

Y/»H J 1* 1^5. *AF ( L ) -154. *A. i ( L 2 - 

1 +61.*AF(L-8)-lQ*^AF (L-10) ) 

: 'd?F ( L )= ( HK^Y2-Y1 ) /HK1 

V 'RtTURT . ■ 


X gM AF ASSEMBLY FLICK 


MO MFSSAUta 


pn f? A ft 0 V F ASSEMBLY 






73 


S T T l : 2 T f ! T 


FORTRAN! SOURCE LIST 


0 

t IS FT 

C K UT T o 

1 


S!60<ruTI0E KUTTA (HA, P”. 

2 


I.C5IC0L LAG 

3 


f I ' : LN S I COAT ( 2 1.0 ) ,$(210. 



101 S (2 10 ) , 02 s ( .210 ) ,A\/( :r 



■- * DK(6),CU 

/, 

A 

L p 

f* U R A T ( 4 1 1 5 . 5 ) 

CALL FLU'iOOCn) 

A, 


CAL 1. FLICK (H,NJ a » AF,DiF,l 

7 


C ALL T ICK { H» UA , S » 0 IS ) 

10 


nci] L = 1 , N] A 

11 


F 3 F { L } =- ( 4F ( L ) *D2 F ( L ) + [ 

3 2 

i I 

02 c ( U =-PK*{ A F ( L ) *D1 S U 

14 


pk ( n=o. 

I r ' 


A '•/ ( 11=1 . 

1 6 


AVV( 1 )=1. 

17 


A V ( 2 1 = 1 . 2 

it!’ 1 - 


A V V (?) = !.. 

21 


A V ( 3 ) = 1 • 

i ? 


A V V (31 = 1.3 

2 3 


A V ! 4 1=1.2 

24 


A V V (4) = 1.3 

25 

2 6 


A V ( 5 1 = I . 5 

AW (5 1 = 1.2 

2 7 


A V ( 6 ) = I . F 

3 0 


AW (A 1 = 1. 2 

3 1 


A V ( 7 1 =1 . 6 

32 


A V 7 { ? 1 = 1 . 6 

3 3 


A V ( 3 ) - 2 » 

34 


AW <? 1=2. 1 

35 


.N$*3 

36 


OL ( 1 1 =0. 

37 


HF ( i. } =(j » 

4 0 


ON (11=0. 

41 


DO { 1 ) =0 . 

47 ' 


I I = I 

14.3 


NK-1 ■ 

4 4 


KK = 1 , 

4 5 


K* 1 . 

46 


feW = l.. 

47 


1 = 1 



C { 1 1 - 0 • 

51 


T ( l)*Cv 

52 


R=Vf,0 

5 3. 

jr / 


P-0. 

* + 

3 3 

.j} ~ fa 

L l) 44 4 1 = 4,5 

56 


CALL PLUM (100) • 

57 


6=0. 

6 0 

6 3 


IF U . 0 r . 2 ) A = l . / 2 . 

I F ( I .ST. 4) A* 1. 

1 6 


Al=nif : (LUA*D2F(U*H 

6 7 


A ? = 1)2 F (1. 5+A*03F(L)*H 


> r » I 1,1 


A 3 = 0 1 S ( L ) + A 4 Q 2 $ ( L ) $ M 
A A = “ F (L } + niF U) *H*A 



FORTRAN SOURCE LIST KUTTA 


| |’':C c: STATS m u j 


7 ? A 0 = 3 (1. 5+01 S(U*H*A 

n ) + l?K { I - 1 . ) 4 4 

74 Y = + A L ( I - I ) * A 

73 ?. = S I - I ) * A 

7 6 I.' •= T f I. )+0\!( I — I }*A 

n v=P+f'c ( i- 1 ) *a 

ICi run ) =H*Y 

1.01 0 L ( I ! =H* Z 

IP? C--‘ C I 5 =-!■'*{ X*A?+ft4*7+B* (U-?.*A1*Y J + ( A5-AI**? ) ) 

i 0 3 r': ( j ) =>-)*</ 

1 06 nill! ) =-H*PR* ( 44*V+A3*X ) 

] 03 4 44 Crj-MINUr 

’ r.7 S ( L + l ) = C ( L ) + 1 . /6. # ( l)K (2)+?.*DK ( 3 ) +2 .* OK ! 4 5 +DKC 5) ) 

j 1.0 P = P + 1./6.*{DL{2)+?.*DL{3) + 2.*!.LC4>+0L(5) ) 

1 i 1 f; = +l./6.*(CM(2)+2.*P v (3)+2. *0 W ( 4) + DM ( 3 ) ) 

si? T ( l + 1 ) = T { L) + l ,/6.*(D^(2)+?.*nr ; . ( 35+2.# 00 ( 4 > + DM ( 5 ) ) 

T 1 3 ft a R + 1 . / 6 . * ( DO {?)+?• *00 ( 35 + 2 .*00 ( 4 ) +00 (5) ) 

\]4 T L = T { L + 1 5 


i IS 
5 16 

J I 7 11 

\ 20 

121 72 

1.22 
I 2 2 
1 . ?6 
L 2 7 
1 30 

131 36 

1.33 
1 3 3 
1, 34 

130 3? 
1 4 0 

141 

142 . 

1,43,. . ■ 

1 4 4 ' . . 

1.4.3 
146 
1.47 
1 30 

131 
1 52 
\ 5 3 
] 54 

} 5 5 ? 5 

1 56 
04 


Vi * J ~ '•:> ! L + 1 ) 

IF (L-NA )31»32,3? 

L — L + 1 

GOTO 3 3 

f f •- j ( •" j 4 ) 

LAG=ABS { P ) . LT. .00001. AND. TT . LT . .00001 

If { LAG ) GO TO . 3 5 

AT ( II ) = TT 

ap( i n*p : 

If (U-NS)36t37»37 
11=1+11 
w*Avtm 

W W * A V V ! If ) 

G0TH3 f • ' ■ 

If (KK .30. 2) 00*035 
r, 1 * ( AT (?) -AT M) ) / t A V i 2 5 - AV ( i) ) 
.C2*(AP(2)-AP(1) )/ t A V ( 2 ) - A V ( 1 ) ) 

C 3 = ( A T { 3 ) -AT (1 5 ) / ( AV V ( 3 ) - AV V U. ) 5 
C4= { AP ( 3 ) *-AP{ 1 ) J / ( AVV { 3 ) — AVV ( 1 ) > 

Of f»=-AP ( 1) 

nr t=- at < 1 ) 

: XC=C 1*C4-C2*C3 

'D A V = f X 4 *0F T- 0 F P *C 3 ) / X C 

I) A V V = f C 1 * OF P -C ? * 0 F T )/XC 


Wb'AV ( 1 ) + DAV 


W W = A V V. ( D+DAVV 


KK=fcK+l 


007 03 3 


RETURN • 


For 

IBMAP ASSEMBLY KUTTA 


0 FfSSAGl 


roi 


a i nvf assembly 



5 


*■’ t 0 


' J : \j 


p i? n 


cno4 
I SM 


Sf.'L'R C e: STATEMENT 


0 HBFTCMm 

1 
7 


3 

4 

5 

6 
7 

10 

12 

1 ? 
1 5 
16 
3 . 7 


1 SUB ( Y*M,H, AREA) 

DIMENSION V ( 2 1 0 ) 

S6: , ~ V = 0. ' ■ 

c i , f; ■ r a 


K = M - i 

004 J = 2 » N, 2 

nc-iI = 3 f K,2 

s u v C c = s u r; o + y ( i > 

* ^ A = * >/ 3. * ( Y { 1. ) +4. * $ !;• 0 V + 2 


RET U 2 
f" ( ; 0 


} x 


FORTRAN SOURCE LIST 


» *SUMGD+Y ( K ) > 

I SNAP ASSEMBLY UB 


N'C; BFSSAOCC FOR ABOVE ASSEMBLY 
04 

is BEING FMTGkFD INTO STORAGE. 


I P-LOR — JOB 000000 



7*r 




APPENDIX 3 



^7 


EGOO 4 


scum 


STATEMENT 


FORTRAN SOURCE 


DIMENSION S T ( 1 5 ) , FN ( 1 0 ) » JC I ( 25 ) ,JCF(25) , JQ( 25! , C{ 20 ) ? 

I AF ( 160,25 ), Q{ 160,25 ), AETA ( 169) ,JE3( 25 ) , 

, (6 J ’ NKV f 10) , SST{ 15,6 ) ,FrN( 19,6 ) ,JJCI (25,6) , JJCF{?5,6) , 

3JJK.H 25 , 6 ) , LCH ( 6 ) ,NI B(6 ) ,NFB (6 > 

LOGICAL LG 

FORMAT (-5 X, 39HS0LUTI0N OGTA I NFD WITH PARA . D I FF . METHOD ) 

FORMAT (2X, 6HAF{ 3 )= , 513. A- , 6H A F(6)=,G14.4> 

FORMAT ( 5 X , 45HS0LUT I ON OBTAINED DIRECTLY FROM ORIGINAL ECO.) 
FORMAT (5X,1 OHPAR AMET>R= , FI 1 .4 ) 


7 

100 

F t.) R M A T ( F 4 . 2 » F 4 . 2 , 2 1 3 ) 

10 

2 09 

FORMAT (11) 

11 

9U4 

FORMAT! UEl 1.4) 

12 

2 3 

FORMAT ( 4HKMP= ? I 3 ) 

13 

14 

FORMAT (2F 5.3) 

14 

15 

FORMAT (71 1 ) 

15 

lb 

FORMAT ( 31 1 ) 

16 

13 

FORMAT ( 3F5. 3) 

17 

11 

FORMAT (3 12 5 

20 

10 

FORMAT ( 2 1 .3 , F5 • 3 5 

21 

117 

FORMA F(F5.3,I3) 

22 


RGA01 17, DB, JA 

24 


R G A 0 1 0 » NK MP * NR » H 

27 


RFAD1 1 , ( NNV ( I ) , 1= 1 ,NKMP ) 

34 


REAOl 1, (NIB< I ) , I=1,NKMP) 

41 


READ1 1, ( MFB ( I ) , I=l,NKMP> 

46 


DO 121 = 1 , NKM’P 

47 


M-NIP ( I ) 

50 


MN = NFL’{ I ) 

51 


M*NNV ( I 5 

52 


READ! 3, ( SST ( K , I ),K=1,M) 

57 


READ14, {FFN!K, I) »K*1,MN) 

64 


READ 15, (JJCI (K, I ) ,K=1 ,N) 

71 


RE ADI 5, ( JJCFtK, I ) ,K = 1 ,N> 

76 


READ15, ( JJEQLK ,1 ) , K = 1 » N ) 

103 : 

12 . 

Continue 

105 


RE AO i 6 , { LCH (K) »K-1, NKMP) 

112 


READ209,M I 

1 14 


AET A { 1 ) =0 . 

115 


NAA=NR-1 

L 16 


D0444K= 1 » NAA 

l 17 

/.f 4 4 

AETA(K+1)=AETA(K)+H 

121 


D017KMP=1 ,NKMP 

122 


LG= ( (LCH(KMP) .E0.2) .OR. (L 

123 


NV=NNV( KMP) 

124 


IOC-NIB(KHP) 

125 


I FE=NF6 ( KMP 5 

126 


001 B J-l , NV 

127 


JCI ( J )*JJCI ( J,KMP ) 

1 30 


JLU( J) = JJf-Q{ J,KMP) 

131 

IT 

JCF( J > = J J C F ( J,KMP) 

l 3 3 


0019 J-i , IBB 

134 

10 

ST ( J)*SST( J, KMP) 

136 


D020C J*l , I F 8 

1.37 

200 

F ; ' ( J ) =F F N( J,KMP) 


7 ^ 


GO 04 
ISN 

141 

14? 

143 

146 

151 

154 

157 

162 

165 

170 

171 

172 
1 73 
1 74 
175 


201 
? 0 2 
20 3 
2 04 
205 
214 
2? 3 

2? *5 
2 26 
2 30 
2 3 1 
004 


SOURCE STATEMENT 


FORTRAN SOURCE LIST 


CALL ANT( NVt JEQ * JO * IC ,MI , JC , JCF ) 

L L M= LCH ( K HP ) 

I F { LCH ( KMP } . E(). 0 ) CALLBST { AF , NR , H } 

uc r !m™;^;c!“aGr ST ' JA,DB ’ H ’ B ’ NV ’ f,R ’ ST ’" , ’ JCF ’ JEU ' JC, ’ J «' 

2KHp B Af!o"c I , jc”Ls’LLM|‘ HI X ,H ’ B ’ NV ’ NR ' S T > >=N, JCF , JEO , JC I , JQ, i C , H I , 

I F { LCH { KMP } . EQ, 3 ) GOT020 
I F ( LCH { KMP ) . EQ . 0 ) GOTO 17 
I F { LCH{ KMP ) . EQ • 1 5GOT021 
I F ( LCH{ KMP) . EQ.2 ) GGT02? 

22 PR I NT 302 

PR I NT23 , KMP 
GO TO 3 04 

2C 0024NS=1,NR 
1 1* t 4 fv vj — 1 , 0 \J 

0 ( MS » N J ) =AF ( AS, MJ ) 

1 Ch5 ;,r:v ri!r:,ST » F ‘‘ ! ’ JCf » JEQ, JCI ,JQ, IC,MI , K M P , A F ,Q,C, JC 

GO TO 3 04 
PRINT 301 
Pk I M T 3 (j 0 f B 
n 0805 N = ]. , N V 

fF( rr/B’rf, ,, PRINT904, ( AF ( I ,NJ , AETA (I ) , 1 = 1 ,nr ) . 

CONTINUE * ? PRINT 904, ( AF ( I , N ) ,AETA ( I ) , 1=1 ,NR) 

GO TO 17 
CONTINUE 
STOP 
ENO 

IBMAP ASSEMBLY 


24 

21 

304 

B 0 5 

t 7 


NO MESSAGES FOR ABOVE ASSEMBLY 


79 


004 

ISM 


FORTRAN SOURCE LIST 


SOURCE STATEMENT 


0 SI8FTC MIST 

1 Sl’P^ri'TINEMIST { JA,DB,H,e,NV,NR,ST,FN,JCF» JEQ, JCI , JQ, IC » Ml , KMP » AF 
1 ,C,C» JC,LG) 

2 LOGICAL LG 

3 DI v CIS SI ON JCF ( 25 } ,JEQ(25),JCI ( 2 5 ) , JQ < 25 ) , JC ( 6 ) , AF { 160 , 25 ) , Q ( 1 

260 » 25), 01( 160,25) »AFF (160) , DIF ( 160) ,C (20) , AF 1 ( 160, 25 ) , ST ( 25 ) , FN ( 

310) 

4 00909 5 =1 » IC 

5 909 C { I )=1 .0000 

7 f''VV = 'V-l 1 . 


10 


IX 53? 1 = 1, MW 

11 


T F { JO ( I ) . CQ« 1 ) G0TD512 

14 


DO 5 1 3M=1 » HR 

1 b 

513 

AFF (N )=AF (N» I ) 

17 


CALL 'TICK (H, NR, AFF, DIF) 

20 


DO 5 1 4 N= 1 , MR. 

21 

514 

AF CM, 1+1 )=01F(N) ' 

2 3 

512 

CONTINUE 

25 


X Y 7 = 0 R / 6 • 

26 


D D 8 = Ij B / 2 • 

27 


b = 0 « 

30 


JJ = 1 

31 

BO? 

DO 80 11=3, NR 

32 


00801 N= 1 » NV 

3 3 

801 

AF l ( I » N ) = A F ( I ,N) 

86 


D=B 

37 


A A = . 1 

40 


CALL RKM { H , D , MV, NR , ST , FN , JCF , JEQ , JC I , JO, I C 
1 , LG 3 

41 


DO 802 N® 1 * NV 

42 


DO 502 1=1, MR 

43 


CALL FLUN (3000) 

4 4 


AF 1 ( I , M ) = A F { I ,N)+Q( I ,N)*UBB 

45 

802 

CONTINUE 

•50 


D=R+D8B 

51 


CALL RKM { H » D » N V * NR , S T , FN , JCF , JEQ » JC I , JQ, IC 
1 , LG ) 

52 


D0803N=1 * NV 

53 


00803 1=1, MR 

: 54: ; 


CALL FLUN (3000) 

55 


AF1 ( I ,N)=AFltI , N) +Q1 ( I , N )*DBB 

56 


on ,M)=0( I ,N) +4.*Q1 ( I ,N) 

57 

803 

CONTINUE 

62 


D=B+DO 

63 


CALL RKM { H , 0,NV , NR, ST , FN , JCF , J EQ , JC I , JQ, IC 
1 , LG ) 

64 


DOB 04 N = 1 ,NV 

6 5 ■ 


00804 1=1 , MR 

66 


QU,N)=Qn,N)+Ql( I,N) 

67 


AF(I,N)=AF(I,N)+XYZ*Q( I,N) 

70 

804 

CONTINUE . '■ ' : 

73 

904 

FORMAT ( 6E15 • 5 ) 

74 


B=8 + 08 ' - • 7 . 

75 


JJ=J J+l 

76 


I F { J J »EQ« JA ) GOTO 808' 


FORTRAN SOURCE LIST M 1ST 


004 

ISN 

101 

102 808 
103 
004 


NO MESSAGES 


SOURCE STATEMENT 

GOTO 8 07 
RETURN 
END 

IBMAP ASSEMBLY MIST 


FOR ABOVE ASSEMBLY. 




%f 

;oo4 


FORTRAN SOURCE LIST 

ISN 


SOURCE STMEfFMT 

0 

$ IB FTC HIX 

1 


SUBROUTINE H I X C H »D , NV , NR , ST, FN, JCF , J EQ , JC I , JQ , i C , MI , KMP , AF, Q , C , J 



1C, LG * LM ) 

2 


LOGICAL LG 

3 


r: V F N S I ON S T { 2 6 ) , JCF { 2 5 ) , JEQ 1 25 ) , JC II 2 5 ) , JO ( 25 ) , JC ( 6 ) , AF 1 1 60 » 2 5 ) , 



10( 160,25) ,FN( 10) ,C(20 ) 

4 


1 F { LM.ECU 2 ) G0T03 

7 


G0T04 

10 

3 

NN = 1 

11 


D01 1=3, NV 

12 


I F { JC I ( I ) • FG* 1 ) G0T02 

15 


GOTO) 

16 

2 

C ( NN ) =AF ( 1 » I 5 

17 


MN = Nf)+l 

20 

1 

CONTINUE 

22 

4 

A A = . 0 1 

23 


CALL RKH H,8,NV»NR,STtFM, JCF, JEQ,JCI ,JQ, IC,M , X M P , 0 , A F » C , J € , A A 


' 

1 , LG ) 

24 


RETURN 

25 


END 

004 


IBMAP ASSEMBLY HIX 

NO MES 

.SAGE! 

i FOR ABOVE ASSEMBLY 


FORTRAN SOURCE LIST 




GOOA 

■ I SN 


SOURCE S' 


ATE v iNT 


0 SIBFTC ANT 


1 


SUBROUTINE ANT ( MV, JEQ » JQ, IC,M I » JC, JCF ) 

2 


DI v Kf:inUJE0{25) » JQ ( 2 5 ) , JC { 6 } , JC F { 2 5 3 

3 


o 

o 

oo 

11 

~z. 

4 

8 

Jc( n=JFwt i') 

6 


- r = ' V- 1 

7 


no 7 i =i f vs 

10 


n= i+r 

1 1 


IF C JGQ ( I) .£0.1) JEC(N)=2 

1A 

7 

CONTINUE 

1 6 


DO 1 1 7N= 1,5 

17 

13.7 

JC ( N )=0 

?l 


IF (MI .EU.0JG0T0116 

2 A 


D01 3 N = 1 , H I 

?:> 

11 

JC { N ) = i 

21 

1 1 6 

I A=0 

30 


DO 91 = 1,NV 

31 


IF ( JCF ( I ) .EG. 0)00 TO 9 

3 A 


l A= I A + 1 

35 

9 

CONTINUE 

37 


I C = I A 

AO 


RETURN 

41 


END 

:ooo4 


I BMAP ASSEMBLY ANT 


M(1 MESSAGES for above assembly 


t 


FORTRAN SOURCE LIST 


% 3 . 


G004 

ISN 


SOURCE STATEMENT 


0 MB FTC RKM 


1 


SUBROUT riFRKr • ! H » R »NV , NR f S* FN, JCF, JEO, JC I , JQ > I C , M I » KMP , U, F , C , JC 


1 

,LC) 

2 


LOGICAL LG 

3 ■ 


LOGICAL LAG, LAD? 6 ) , , , „ -_ 4 

4 


n l*!l:flSI QMS ( 2 >) , FN C 10 ) ,JCF!25) , JEQC25) »JCI ( 25 ) » JO ( 25 ) ? F C 160,25) 


7 

X 

U ( 1 6 0 » 2 5 ) » 3Q ( 5 ) » 0 ( 1 1 » 1 1 ) » C (■ 20 ) » G 1 1 1 ) t DF { 5 » 2 5 ) * F F { 1 60 ) »DFF ! 160) 


u- 

Yiil , 11] ,1(11) ,LA(11, 11 ) »CY! 11 til) tXdl) » JC( 6 ) ,LSUM! 11 ) 

5 

200 

FORMAT (£14.4) 

6 

212 

FORMAT (2X,4HG! I > ,£14.4) 

7 

?11 

FORMAT (413) 

1 0 

213 

FORMAT ( ?X » I 3 , 7HF ( 1 * 1 1 ) , 2 1 3 ) 

11 

214 

FORMAT (2X,2HDF, £14.4,213) 

12 

21 5 

FORMAT ( 2X , 8 HY { J , NN ) = » E14 .4 , 41 3 > 

13 

216 

FORMAT (2 X , 9HLA { IP, IN }»E 14*4, 213) 

]> 

217 

FORMAT (2X »8HLSUM( IP) ,213) 

15 

212 

FORMAT (2X, 8 HCY ( M » LL ) ,E 14.4,2 13! 

16 

219 

FORMAT ( 2X,8HCY( I , JK ) , £14. 4, 213) 

17 

220 

FORMAT (2X,4HX II ) ,2E14.4) 

20 

2 21 

FORMAT ( 2X ,4HC ( I ) » 2E14.4 ) 

21 

222 

FORMAT! 2X , 2H I = , 13) 

22 

400 

FORM AT { 1 X , 2HDF , 7E14.4 ) 

23 

33 

FORMAT ( 2F5 . 2 ) 

24 

32 

FORMAT! 3F5. 2) 

25 

31 

FORMAT! 1 0 X , 7 1 2 ) 

2 fo 


CALLFLUN! 3000) 

27 


PR I NT 31 , ( JEO! I ) , 1 = 1 »NV) 

34 ■ 


PRINT31,! JCF (I), 1=1, NV) 

41 


PRINT 31 , ( JCI ( I ) » 1=1 »NV) 

46 


PRINT31 , ! JO! I), 1=1, NV) 

53 


PR INT32 , S ( 1 ) , S ( 2 ! , S ( 3 ) 

54 


PR I NT 33, FN! 1 ) ,FN!2) 

55 


PR I NT 200 , A 

56 


PR I N T 2 1 1 , KMP , N V » I C .» M I 

57 


m m =. 1 

60 


By 12) *.001 

61 


BO ! 3 ) = * 0001 

62 


BO! 41*. 00001 

63 


8 Q( 5 ) =,000001 

64 


PRINT221, C ! 1 ) , C 1 2 ) 

65 


ISUM=0 

66 

461 

Z=A*2. 

67 


KK = 1 

70 


KKK = 1 

71 


IF ( IC.EG.0)G0TC418 

74 


0010 1=1* IC 

75 


D010 J=lt IC 

76 


IFt I.EO.-J)D(ItJ)“Z 

101 


IF! l.NE»J)DC I , J ) = 0 . 

104 

10 

. CONTINUE , 

1 C 7 

418 

JK= IC+i 

110 


KK = l 

ill 


KKK = 1 

1 12 


ZU=1 .00000/6.00000 

113 

100 

00320 J=l, JK 


FORTRAN SOURCE LIST RKM 




GOOA 


ISN 


114 


115 


120 


121 


124 


127 

11 

131 


132 

4 0 

133 


134 


135 


14 0 


143 


146 


147 


150 

13 

151 


15? 


1 5 3 

15 

1 54 

12 

156 


157 


160 

19 

162 


16 3 

547 

165 

21 

166 


167 


3.70 


171 


174 


177 


202 


205 


210 


213 


2lh, 


217 

501 

22:0 


221- 


222 

20 


225 

22b 

22.1 


232. 

2 35 739 

237 

240 

243 

244 300 
2 45 

250 

251 
2 54 


SOURCE STATEMENT 


CMLFl.UU(30bO) 

IF ( I C • EC . 0 ) GOT 04 C 

0011 1 = 1, IC 

I F { J « EO. 1 ) 0 ( I ! =C { I ) 

!HJ.GT,l)G!I)=Cm+D(J-l,I) 

CONTI Xi lJF 

PSINT221 ,G( 1 ) ,G(2) 


i • •' .1 

0012 I Z=1 ,NV 
I F { JOIIZ) .EC.DGGTC15 
I F { JC I { I Z ) . EQ.O ) GQTOl 3 
I F C JCH IZ ).EG.1)F(1,IZ)=G(N8) 

N R = N R + 1 
GOTO 1 2 

FI 1 , IZ)=S(NS5 
NS=NS+1 
GO TO 12 
F ( 1 , 1 L ) =0 . 

CONTINUE 


K= 1 

0019 M= 1 , NV 
DF ( 1 » M ) =0 
DO 54 7 M= 1,5 
OF ( M , 1 ) =0 . 

DO 20 M=2,5 
M0 = 0 

DO 20 MA=2,NV 
CALLFLUNt 500) 
IF(M.EG.2)A=0. 


IF ( ?'• « FC * 3 ) A =(' . 5 CO 00 
IF (H.EQ.4)A=. 50000 

IF ( H. EQ» 5 1 A= 1 .00000 . 

IF ( JFC( HA ) . EQ.O IDF ( M, MA) = H*(DF { M-l ,MA+1 )#A+F{K,MA) ) 

I F ( JEC ( MA ) * EQ • 1 ) GOT0501 

IF CJFOI MA ) . EQ.2 5 DF { F ,MA ) =DF ( 1 » MA ) 


GO TO 20 , 

MO=NO+ 1 

CALL FUNSON(H,B,A,K,NV, JO, H, MA ,Q, F ,DF,DFP, KMP ,N0 ) 
DFtM,«A)=DFP 
CONTINUE 
D0739M=1 , NV 

lF(JO^MK£e!o)F{K+l,M)=F(K,M)+ZU*(DF(2,M+l)+2.«DF(3,M+l) 

1 +2»*DF( 4, M+l ) +DF ( 5, M+l ) 3 

I F t J C ( M ) . EO » 1. ) F t K+ 1 , M > =0 • 


CONTINUE 

4» j[ 

IF { K. 6Q.NR I GOT0300 


GO TO 2 l 
NN*1 

IFUC.E0.0)GOTO410 


0023 M= 1 , N V 

IF ( JCF(M) ,E0.0)G0TQ23 
I F ( JCF( M) .EO. 1 ) Y ( J , NN )=F ( NR, M) 



1 iEG004 


I $N 


SOURCE 


TCMEMT 


SOURCE LIST RKH 


257 
2 60 
2 61 

263 

264 
2 65 
2 66 

2 67 
272 
275 

3 CO 
301 
.302 
3 03 

304 

305 

310 

311 
313 
316 
320 
3? I 
324 
327 
332 
335 
340 
34 3 
346 

351 

352 
3 53 

354 

355 

356 

36 1 

362 

363 

364 

366 

367 

3 74 

401 

402 

4 C 3 
405 

412 

413 

414 

415 
' 416 

421 
4 22 
42 3 
423 


PR I NT 2 1 5 » Y ( J , MN ) ,J,!*N,NR t M 

3;>» — -'j T,’ 1 

2 3 CONTINUE 

704 ;U I P= 1 » 5 

304 01 I N=i » IC » 

T( I “j } =ABS ( Y ( J , I N ) -FN I IN) ) 

L4G = T (IN) »L T .8Q( IP) 

IF (LAG) LA( IP,IN) = i 
IF? .NOT . LAG ) L A ( IP, IN)=0 

4 u i CONTINUE 

DO 40 2 I P= 1 , 5 
I SUM=0 

00402 IU=1, IC 
ISUM=LA{ IP, IN)+ISUM 
LSUM? IP ) = I SUP 

402 CONTINUE 
D0403 I P= 1. , 5 

403 LA!) { !P)=LSUM( IP) .EO. IC 
I F •( L A U ( 5 ) ) G0TC410 

320 CONTINUE 

LAG=KK.FQ.2 

IF? ( . NOT . LAD ( 1 ) ) . AND. (LAG). AND. { JC ( 1 ) .EG. 1 ) )GDT0460 
IF? ( . NOT. LAD ( 2 ) ) .AND. (LAG). AND. ( JC { 2 ) .EO. 1 ) )G0T0460 
IF ( ( . NOT. LAD (35). AND. ? LAG) . AMD. ( JC { 3 ) .EQ. 1 } ) G0T0467 
IF ( (.NOT. LAD (4) ) .AND. ( L AG ) . AND . ( JC ( 4 ) .EQ.l) )G0TO462 
IF (LAD? 5) .AND. ( JC (5 ) .EO. 1 )) GOT 04 10 

IF C ( . NOT . LAD ( 5 ) ) .AND. { L AG ) . AMD . ( JC { 5 ) .FQ.l) )G0T0463 

If (LAD (4) .AMD. < JC(4) .EQ.l) )GOT0410 

I F ( K K . 6 Q . 2 ) G0T04 1 0 

DO 2 5 L=2, JK 

D025P=1, IC 

LL = L- 1 

CY(M,LL)=(Y(L,P)-Y( 1,M ) Ml 
PRINT218,CY(M,LL) » M , LL 

25 CONTINUE 

DO 2 6 1 = 1, IC 
CVM,JK) = FN(n-Y(l,n 
PR I N T 2 1 9 » € Y ( I , J K ) , I , J K 

26 CONTINUE 

CALL ALEOt JK,CY,X) 

PR I NT 220 , (XU ) ,1 = 1, IC ) 

PR I NT 22 1 » CC1I) ,1 = 1, IC) 

D07QI =1 » I C 
C( I )=C( I)+X( I ) 

70 CONTINUE 

PRIMT221, (C(I),I=1,IC) 

'kk=kk+i 
JK = 1 

GOTO ICO " 

410 0041=2, NV 

IF( JCm.EO.l) GOT05 
G0T04 

5 D06N=1,NR 

6 FF (N)»FtN, 1-1 ) 

CALL TICK (H,NR,FF, OFF ) 



EG004 

] ISM 


SOURCE STATEMENT 

1 426 


D07N=1,NR 

| 427 


FIN, I ) = DFF ( N ) 

I 4 30 

7 

CONTINUE 

: l 4 3? 

4 

CONTINUE 

| 4 34 


GOT 04 79 

| 4 35 

460 

I F ( L G ) A = . Q 0 5 

1 440 


A = 0 . 0 2 5 

441 


GOTOA61 

442 

4 62 

I F I LG ) A=* 002 

4 4 5 


A= . 005 

- ‘t 4 6 


G0T0461 

447 

46 3 

IF { LG ) A = « 001 

; 4 6? 


A= • 00? 

45 3 


G0T0461 

| 4 54 

479 

RETURN 

i ; f 45 5 


END 

|EG004 

i 

MO MES 

SAGES 

FOR ABOVE ASSEMBLY 


24 

FORTRAN SOURCE LIST RK.f" 


IBMAP ASSEMBLY RKM 


\ 





F b " Ui+ FORTRAN SOURCE LIST 

I I 5N SOURCE STATEMENT 

; I 

I 0 ilBFTC TICK 

i 1 SUBROUTINE T I CK ( H , NA , AF » DI F ) 

2 DIMENSION AFI210) , DIF (210) 

3* CALLFLUNI 3000 5 

4 NNA=NA-4 

5 NB=NNA+ I 

6 HH=1 ./( I2.*H) . 

7 HJ=l./( 12.*2.*H) 

10 HK={ 1./2. )**4 

11 HK I=HK- I . 

12 D046L=I ,4 

13 X1 = I-25.*AF(D+48.*AF CL+1 }-36.*AF(L+2 ) +16.*AF ( L+3 1-3. *AF ( L *4) )*HH 

14 ■X?*HJ*C-2'5.*AF(U+48.*AFlL , + 2>-36.*AF(L+4}+16.*AF< L+6 ) -3 . * AF ( L+G ) ) 

15 46 OIFIl )=(HK*X2-X1)/HK1 

17 D047L=5,MMA 

20 X1=HH*( AF (L-2)-B**AF(L-i)+6.*AF(L+l )-l.*AF(L+2) ) 

. 21 X2=HJ*{AF(L-4>-8.*AF(L-2)+8.*AF(L+2}-l.*AF(L+4)) 

22 47 D1F(L)=(HK*X2-X1)/HK1 

24 D046L=NP f NA 

25 X1=HH*(>5.*AF( L)-48.*AF (L-l ) +36.*AF ( L-2) -16 .*AF ( L-3 ) +3. *AF l L-4) 1 

?! 26 X2=HJ*(25.*AF<L)-48.*AFlL-2)+36.*AF{L-4)-16.*AF(L-6)+3.*AF(L-8> ) 

| 27 48 OIF { L ) = (HK*X2-X1 ) /HK1 

I 31 RETURN 

i 3? END 

SEGO04 I BMAP ASSEMBLY TICK 


j 00 MESSAGES FOR ABOVE ASSEMBLY 






1 NO MESSAGES FOR AOOVfi ASSEMBLY 

I 


';tff=G 0 D 4 

I SN 


SOURCE STATEMENT 


FORTRAN SOURCF LI 


N'EGOf 


$ IBFTC B ST 

SUB R OUT I ME B ST ( A F , NR t H ) 

DIMENSION AF( 160,25) 

E T A = 0 . 

002211=1, NR 

AF( I »l)=ETA+EXP(-ETA)-l. 

AFI I ,5) =1. + (1.5-1. ) *EXP { -ET A ) 

221 ETA = f: TA+H 
RETURN 
END 

I8MAP ASSEMBLY 


NO MESSAGES FOR ABOVE ASSEMBLY 

MFC DO 4 


t BLUR — JOB 


MEMORY 


) INCLUDING IOCS 

JOCK ORIGIN 

\ 

jlBBER OF FILES - 

1 S.FBIN 
\ S.FPOU 

i PROGRAM 


12220 

12243 


■] 

l 

DECK 

f 



12266 

l 

DECK 

‘ MIST 



35051 

K ■ 

DECK 

'MIX 



56135 

J 

■k& c 

DECK 

' ANT 



56342 

i 

DECK 

1 RKM " 



56546 

I 

DECK 

'* MEG- 



63161 

1, 

deck 

’TICK 



6 3704 

i. 

DECK 

•FUN SON 



64371 

i.' 

DECK 

’ BST : 

f • • 


65025 

u : 

SURR 

' INSYFB 



65140 

i. 

SUBR 

•0USYF8 


: [ 

65177 

*. 

SUBR 

•POSTX 



65230 

u 

i. i. 

SUBR 

•CNSTNT 



65542 

SUBR 

• FPR 



65551 

k. 

SUBR 

» FRO 



65552 

'f6# 

SUBR 

•ICS 



65553 

■ fc- 

SUBR 

* RWD 

*' - 


66032 


SUBR 

»f:CV 

t 


67206 


SUBR 

* F C V 



674 54 

SUBR 

• HCV 



67546 

ti ■ 

SUBR 

•ICV 



67651 

■ & 

SUBR 

• xcv 

1 ' 


'■ 67671 

?3, 

SUBR 

* I. NT J 



. 67707 

: I#.- 

■ SUBR 

* FFC 



70223 

mm 

SUBR 

* FPT 



70645 


00000 THRU 12211 
12220 


12266 THRU TIT 13 



